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ABSTRACT 

We solve the one-dimensional neutrino-heated non-relativistic magnetohydrodynamic (MHD) wind prob- 
lem for conditions that range from slowly rotating (spin period P > 10 ms) protoneutron stars (PNSs) with 
surface field strengths typical of radio pulsars (B < 10 13 G), to "proto-magnetars" with B as 10 14 - 10 15 G in 
their hypothesized rapidly rotating initial states ms). We use the relativistic axisymmetric simulations 

of Bucciantini et al. (2006) to map our split-monopole results onto a more physical dipole geometry and to 
estimate the spindown of PNSs when their winds are relativistic. We then quantify the effects of rotation and 
magnetic fields on the mass loss, energy loss, and thermodynamic structure of PNS winds. The latter is partic- 
ularly important in assessing PNS winds as the astrophysical site for the r-process. We describe the evolution 
of PNS winds through the Kelvin-Helmholtz cooling epoch, emphasizing the transition between (1) thermal 
neutrino-driven, (2) non-relativistic magnetically-dominated, and (3) relativistic magnetically-dominated out- 
flows. In the last of these stages, the spindown is enhanced relative to the canonical force-free rate because of 
additional open magnetic flux caused by neutrino-driven mass loss. We find that proto-magnetars with P w 1 
ms and B > 10 15 G drive relativistic winds with luminosities, energies, and Lorentz factors (magnetization 
a ~ 0.1 - 1000) consistent with those required to produce long duration gamma-ray bursts and hyper-energetic 
supernovae (SNe). A significant fraction of the rotational energy may be extracted in only a few seconds, 
sufficiently rapidly to alter the asymptotic energy of the SN remnant, its morphology, and, potentially, its nu- 
cleosynthetic yield. We find that winds from PNSs with somewhat more modest rotation periods (w 2- 10 
ms) and with magnetar-strength fields produce conditions significantly more favorable for the r-process than 
winds from slowly rotating, non-magnetized PNSs. Lastly, we argue that energy and momentum deposition by 
convectively-excited waves may be important in PNS winds. We show that this further increases the likelihood 
of successful r-process, relatively independent of the PNS rotation rate and magnetic field strength. 

Subject headings: stars: neutron — stars: winds, outflows — supernovae: general — gamma rays: bursts — 
stars: magnetic fields — nuclear reactions, nucleosynthesis, abundances 



1. INTRODUCTION 

On timescales < 1 s following the core collapse of a mas- 
sive star, neutrino emission from the resulting hot, delep- 
tonizing protoneutron star (PNS) may play an essential role 
in launching the supernova (SN) shock (Herant et al. 1994; 
Burrows, Hayes, & Fryxell 1995; Janka & Miiller 1995) or 
in generating large-scale anisotropies through hydrodynami- 
cal instabilities (Blondin et al. 2003; Scheck et al. 2006; Bur- 
rows et al. 2006a,b). Independent of how the explosion is 
initiated at early times, a small fraction of the PNS's cool- 
ing neutrino emission continues heating the surface layers 
of the PNS, driving a persistent thermal wind into the cav- 
ity evacuated by the rapidly-expanding SN shock (Duncan, 
Shapiro, & Wasserman 1986; Woosley et al. 1994); this post- 
explosion neutrino-driven mass loss persists for the duration 
of the PNS's Kelvin-Helmholtz cooling epoch, which lasts a 
time tkh ~ 10— 100 s (Burrows & Lattimer 1986; Pons et 
al. 1999). 

Once the SN shock has been launched and the PNS cool- 
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ing epoch has begun, neutrino-driven winds from non-rotating 
non-magnetic PNSs are unlikely to be energetically impor- 
tant on the scale of the kinetic energy of the accompany- 
ing SN (£sn ~ 10 51 ergs). However, in the presence of a 
sufficiently strong global magnetic field the dynamics of a 
PNS's neutrino-heated outflow are significantly altered (e.g., 
Thompson 2003a,b). This point is germane because ~ 10% 
(Kouveliotou et al. 1994; van Paradijs et al. 1995; Lyne et 
al. 1998) of young Galactic neutron stars possess significantly 
stronger surface magnetic fields (~ 10 14 - 10 15 G) than those 
usually inferred from pulsar spin-down estimates ("magne- 
tars"; for a recent review see Woods & Thompson 2004). 
While the precise origin of these large field strengths is un- 
certain, it has been argued that their amplification occurs via 
a dynamo during tkh (Duncan & Thompson 1992; Thompson 
& Duncan 1993). The efficiency of the dynamo is determined 
in part by the core's initial rotation rate Slo = 2tt/Pq, and ana- 
lytic arguments suggest that the formation of global magnetar- 
strength fields might require Po ~ 1 ms rotation at birth. Such 
rapid rotation would also alter the dynamics of the PNS wind 
and provide a reservoir of rotational energy significant on the 
scale of the accompanying SN explosion: 

E mt ~ 2 x 10 52 M lA R 2 w P^ergs, (1) 
where Mi .4 is the PNS mass in units of 1.4 Mq, R 10 is the 
radius of the PNS in units of 10 km, and P ms is the initial PNS 
rotation period in ms. 

Previous authors have suggested that if magnetars are in- 
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deed born rapidly rotating, their rotational energy could be 
efficiently extracted through a magnetized wind (Usov 1992; 
Thompson 1994; Wheeler et al. 2000; Thompson, Chang, & 
Quataert 2004, henceforth TCQ). Indeed, a proto-magnetar 
wind's energetics, timescale, and potential for highly rel- 
ativistic outflow resemble those of the central engine re- 
quired to power long-duration gamma-ray bursts (LGRBs). 
TCQ argue that an accurate description of magnetized PNS 
spin-down must include the effects of neutrino-driven mass 
loss. The significant mass loss accompanying the Kelvin- 
Helmholtz epoch may open the otherwise closed magneto- 
sphere into a "split monopole"-like structure, enhancing the 
early spin-down rate of the PNS. In addition, the mass load- 
ing of a PNS wind, and hence its potential asymptotic Lorentz 
factor T, is largely controlled by the neutrino luminosity dur- 
ing the Kelvin-Helmholtz epoch. This "baryon-loading" is- 
sue is particularly important in the LGRB context, where 
T ~ 10- 1000 are typically inferred. 

Because magnetar births are relatively frequent, most can- 
not produce classical LGRBs, which only occur in ~ 0. 1 - 1 % 
of massive stellar deaths (e.g., Paczynski 2001, Podsiadlowski 
et al. 2004, Piran 2005); however, a more common obser- 
vational signature of magnetar birth may include less ener- 
getic or mildly relativistic events, which could be observable 
as X-ray transients or unusual SNe. Under some circum- 
stances, the asymmetric energy injection from proto-magnetar 
winds could produce global anisotropies in the SN remnant, 
as has been detected through polarization measurements of 
some SNe (e.g., Wang et al. 2001, 2003). In addition, if sig- 
nificant rotational energy can be extracted sufficiently rapidly 
following the launch of the SN shock, the nucleosynthetic 
yield of the SN could be altered (TCQ), which could ex- 
plain nickel-rich, hyper-energetic SNe such as SN 1998bw 
(Galama et al. 1998) or SN 2003dh (Hjorth et al. 2003; Stanek 
et al. 2003). Lastly, LGRBs occurring without associated SNe 
(e.g., GRBs 060505, 060614; Fynbo et al. 2006; Gal- Yam et 
al. 2006; Delia Valle et al. 2006) may be accommodated if 
a proto-magnetar results fro m the accretion-induced collapse 
(AIC) of a white dwarf (see §04). 

Quite apart from the possible impact of PNS winds on 
the surrounding SN shock, PNS winds themselves have of- 
ten been considered a promising site for the production of 
r-process nuclides. However, the conditions necessary for a 
successful third peak (A w 195) r-process have not been real- 
ized in detailed studies of non-rotating, non-magnetized PNS 
winds (e.g., Qian & Woosley 1996; Cardall & Fuller 1997; 
Otsuki et al. 2000; Thompson et al. 2001). Given the, as 
yet, unidentified site of Galactic r-process enrichment and the 
relatively large birthrate of magnetars, it is essential to con- 
sider what effects a strong magnetic field and rapid rotation 
might have on nucleosynthesis in PNS winds. Conversely, 
if the nucleosynthetic yield from rotating, magnetized PNS 
winds could be well-determined theoretically, the birth period 
and magnetic field distribution of neutron stars could perhaps 
be constrained from observations of heavy elemental abun- 
dances. 

1.1. Stages of PNS Evolution 

A sufficiently-magnetized PNS wind goes through at least 
three distinct stages of evolution following the launch of the 
SN shock (f = 0): 

(1) During the earliest phase the PNS surface temperature 



is so high (> 5 MeV) and the radius of the PNS is so large 
(R v ~ 20 -50 km) that, even for a magnetar-strength field, the 
post-explosion outflow is likely to be purely thermally-driven. 

(2) By t ~ 1 s the SN shock has propagated to well out- 
side the sonic point of the PNS wind. The PNS will con- 
tract and cool to a point at which, if the surface field is suffi- 
ciently strong, the outflow becomes magnetically dominated. 
The strong magnetic field enhances angular momentum and 
rotational energy loss by forcing outgoing fluid elements to 
effectively corotate with the PNS surface out to the Alfven 
radius (Ra) at several stellar radii, in analogy with classic 
work on non-relativistic stellar winds and the solar wind (e.g., 
Schatzman 1962, Weber & Davis 1967, and Mestel 1968). For 
sufficiently rapidly rotating PNSs, the spin-down timescale 
(tj = Sl/Sl) can be comparable to tkh, implying that much 
of the rotational energy of the PNS can be extracted in a non- 
relativistic, but magnetically-dominated wind. The neutrino 
luminosity at these relatively early times is still large (e.g., 
L v ~ 10 52 ergs s" 1 ; see Pons et al. 1999, Fig. 14), and, for 
sufficiently rapid rotation, mass loss is significantly enhanced 
by centrifugal flinging (TCQ). The outflow during this phase 
will be collimated about the PNS rotation axis by magnetic 
stresses (e.g., Bucciantini et al. 2006; hereafter B 06). 

(3) As the PNS continues to cool, the luminosity and mass- 
loading decrease to the point at which corotation is sustained 
out to nearly the light cylinder (R^ = 2nc/P ~ 48 P~' km). 
The spin-down rate then becomes approximately independent 
of the mass loss rate M and the flow becomes relativistic, ob- 
taining high magnetization a = <f>|il 2 /Mc 3 , where is the 
total open magnetic flux per 4ir steradian (Michel 1969). Al- 
though an ultra-relativistic, pulsar-like wind is inevitable soon 
after tkh because M abates as the neutrino luminosity van- 
ishes, relativistic outflow can begin prior to the end of the 
Kelvin-Helmholtz phase. Because this outflow is accompa- 
nied by significant mass loss, it is generally only mildly rel- 
ativistic. Although the energy extracted via relativistic winds 
is primarily concentrated at low latitudes (B06), the confining 
pressure of the overlying, exploding, stellar progenitor (e.g., 
Wheeler et. al 2000; Uzdensky & MacFadyen 2006) or the 
walls of the collimated cavity carved by the preceding non- 
relativistic phase may channel the relativistic outflow into a 
bipolar, jet-like structure. 

The primary focus of this paper is to delineate the magnetic 
field strengths and rotation rates required to significantly al- 
ter the characteristics of early PNS evolution by solving the 
one-dimensional (ID) non-relativistic MHD, neutrino-heated 
wind problem for conditions that range from normal pulsars 
to proto-magnetars. In particular, we quantitatively explore 
the transition that occurs between stages (1) and (2) above; 
with these results we analyze some of the immediate conse- 
quences of neutron star birth. We defer a detailed study of the 
transition between stages (2) and (3) to future work, but we 
do explicitly address the parameter space of the a = 1 bound- 
ary. Relativistic, mass-loaded MHD winds have been stud- 
ied recently by B06 in two dimensions for both a monopo- 
lar and aligned-dipolar field structure, assuming an adiabatic 
equation of state. While the work of B06 is critical to under- 
standing the multi-dimensional character of PNS winds (for 
instance, the degree of collimation and the fraction of open 
magnetic flux) it does not address the neutrino microphysics 
necessary for a direct application to PNS environments; this 
work and that of B06 are thus complementary in studying 
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PNS spin-down in the presence of rapid rotation and a large 
magnetic field. 

1.2. This Paper 

In Sj2]we enumerate the equations of MHD ( £12.11 ), discuss 
the relevant microphysics ( 32.21 i, and elaborate on our numer- 
ical methods ( 32.3t . Section [3] presents the results of our 
calculations and examines the regimes of magnetized PNS 
wind evolution. Section |4] discusses the implications of this 
work, examining the time-evolution of a cooling, magnetized 
PNS weighing the implications for LGRBs and hyper- 

energetic SNe ( 34.2| >, and considering the viability for third- 
peak r-process nucleosynthesis in magnetized, rotating PNS 
winds ( 34.31 l, In 34.3. ll we consider the effects that wave heat- 
ing have on the r-process in PNS winds and in 34.4l we briefly 
discuss other contexts in which our calculations may be ap- 
plicable, including the accretion-induced collapse of white 
dwarfs, "collapsars," and merging neutron star binaries. Fi- 
nally, section[5]summarizes the conclusions of our work. 

2. PNS WIND MODEL 

2.1. MHD Equations and Conserved Quantities 

Making the simplifications of Weber & Davis (1967), we 
restrict all physical quantities to be solely functions of radius 
r and time t, and confine our analysis to the equatorial plane 
so that the magnetic field B = [B r ,Z?0] and fluid velocity v = 
[v r , vj,] have no 9 components. We employ Newtonian gravity 
for a PNS of mass M acting on gas of density p. We also as- 
sume that the plasma is a perfect conductor with an isotropic 
thermal pressure P. Under these restrictions the time evolu- 
tion equations of non-relativistic MHD are 



dp 



1 d 



dt 



dr 



dVr_ _ V J_ _ d\r_ 

dt r r dr 
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p dr 



4np 



r dr 
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(5) 



<9v0 1 I" B r d 
dt r 4irp dr 

dB^ _ 1 d 
dt r dr 

We have neglected neutrino radiation pressure in equation ([3]l 
because the neutrino luminosity is always below the neutrino 
Eddington limit. In steady state, equation (f2]i gives a radially- 
conserved mass flux: 



M = AQr 2 pv r , 



(6) 



where Aft is the opening solid angle of the wind. While the 
evolution equations considered are formally valid only in the 
equatorial plane, quoted values for M will be normalized to 
Ail = 4tt, as if the solutions were valid at all latitudes. As 
discussed in 34.11 this normalization is an overestimate be- 
cause the PNS's closed magnetic flux prevents mass outflow 
from near the equator and centrifugal flinging concentrates M 
at relatively low latitudes (B06). 

We take the radial magnetic field structure to be that of a 
"split monopole": B r = B v (R v jr) 2 , where R„ is the radius 
of the PNS and B v is the monopole surface field strength; 



this prescription conserves the magnetic flux 47t$b = 4-Kr 2 B r . 
B06 show that the spin-down in dipole simulations can be ex- 
pressed in terms of an equivalent monopole field, which de- 
pends on the fraction of open magnetic flux. As discussed 
further in 34.11 we can therefore relate our monopole spin- 
down calculations to the more realistic dipole simulations of 
B06. 

In steady state, manipulation of equations © and (0 gives 
the conserved specific angular momentum C and "the conse- 
quence of induction" X (e.g., Lamers & Cassinelli 1999): 



£ — £gas "t" ^mag — Weft ' 

1= r{v 4> B r -v r B 4 



rBrB^ 
4npv r 



(7) 
(8) 



The total rate of angular momentum loss from the PNS is 
therefore J = CM. 

For conditions of interest, photons are trapped and advected 
with the wind, providing no significant energy transport on the 
timescales of interest (Duncan, Shapiro, & Wasserman 1986). 
Instead, the PNS evolution is controlled by its neutrino lumi- 
nosity L v , which provides heating q[, (per unit mass) above 
the PNS surface. Including net neutrino heating (q v = q*—q~; 
see 32.21 i provides a source term in the wind entropy equation: 



dS . 
T lt =q - 



(9) 



where d/dt = d/dt + v r (d/dr), S is the wind entropy per unit 
mass, and T is the wind temperature. The asymptotic wind 
entropy S" quantifies the total net heating a parcel of gas ex- 
periences as it is carried out by the PNS wind. As discussed 
further in 32.21 when applying equation (O we assume that 
the wind's composition (electron fraction) remains constant. 

With net heating the Bernoulli integral B is not constant 
with radius (in steady state); instead it receives an integrated 
contribution from q v : 



MAB = M{B(r)-B(r )}= / pq v AVtr' 2 dr 



where 



B= l -(v 2 + vl) + h-™ + £ m . ig n, 



(10) 



(11) 



h = e + P/p is the specific enthalpy, e is the specific internal 
energy, £ m ag^ is the specific magnetic energy, and il is the 
stellar rotation rate. 

To assess the impact of the PNS wind on its surroundings, 
it is useful to define r\, the ratio of the rotational power lost by 
the PNS, E ml = QJ, to the asymptotic wind power E a = B a M: 



QJ _ nc 

77 = ~E"~ 



(12) 



where B a is the Bernoulli integral evaluated at large radii. 
Magnetized winds are typically assumed to have rj w 1, but 
we find that in some cases 77 ^> 1 because much of the rota- 
tional energy is used to unbind the wind from the PNS (see 



2.2. Microphysics 

The local heating and cooling rates relevant to the Kelvin- 
Helmholtz cooling phase have been extensively evaluated in 
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efforts to quantify unmagnetized, slowly-rotating PNS winds 
as an astrophysical site for r-process nucleosynthesis (e.g., 
Qian & Woosley 1996, Thompson et al. 2001; hereafter TO 1). 
For our calculations we adopt the heating and cooling rates 
used in the analytic work of Qian & Woosley (1996, here- 
after QW) for v jv annihilation and for the charged-current 
processes v e + n <-> p + e~ and v e + p «-> n + e + . For heating 
and cooling from inelastic neutrino-lepton scattering we use 
the rates of T01. The dominant contributions to the net heat- 
ing rate are the charged-current processes, but scattering and 
annihilation become more significant as the entropy of the 
wind increases and the thermal pressure becomes radiation- 
dominated. The scattering and charged-current rates effec- 
tively vanish once T — > 0.5 MeV because the e~/e + pairs an- 
nihilate and the nucleons combine into a particles. We ar- 
tificially take this cutoff into account by setting q u = for 
T < 0.5 MeV. 

Because the charged-current interactions modify the neu- 
tron abundance of the outflow the electron fraction Y e should 
be evolved in addition to the temperature and pressure. How- 
ever, the dynamics of the wind are not sensitive to the precise 
profile of Y e and thus, for simplicity, we take Y e to be fixed at a 
reasonable asymptotic value at all radii: Y e = Y" = 0.4. This is 
a good approximation given how rapidly in radius Y° obtains 
in non-rotating, unmagnetized calculations (see T01, Fig. 7) 
and how relatively weakly the total heating and cooling rates 
depend on Y e for Y e < 0.5. 

As in T01 we calculate local neutrino fluxes by consider- 
ing a single (for all neutrino species) sharp, thermal neutri- 
nosphere at R v . To the desired accuracy of our calculations 
this approximation is generally good, even when the density 
scale height (and hence nucleon-absorption optical depth) is 
extended by rapid rotation. For this paper we index stages of 
the PNS thermal evolution in terms of the anti-electron neu- 
trino luminosity L Pc . We scale all other neutrino luminosi- 
ties (L Vt ,L v ^,L Pli ,L VT , and L 9t ) as in TCQ: L Ue =L Q J\3 = 
1 .08L„ M , where /i denotes each of the other four neutrino/anti- 
neutrino species. Note that the total neutrino luminosity is 
then L v ~ 4.6L Pf . Following T01, all first energy moments at 
the neutrinosphere ((e v ) = (E^)/(E l/ ), where E v is the neu- 
trino energy) were scaled with luminosity as (e v ) oc L l J 4 , an- 
choring { (e„, ) , (e Pe ) , (e„„ ) } at { 1 1 , 14, 23}Me V for L Pe>51 = 8, 
where Lp f ,5i is the anti-electron neutrino luminosity in units 
of 10 51 ergs s" 1 . Higher energy moments necessary for the 
heating calculations ((e" ), (eg ), etc.) are related to the first 
through appropriate integrals over the assumed Fermi-Dirac 
surface distribution. We should note that the relationship be- 
tween the neutrino luminosity and mean neutrino energy we 
have assumed, while a reasonable approximation, is likely to 
be more complicated. For example, Pons et al. (1999) find 
that the mean energy is roughly constant for the first ~ 10 s of 
cooling, despite the fact that the neutrino luminosity decreases 
monotonically (see their Fig. 18). 

The neutrino heating and cooling rates discussed above will 
be modified by the presence of magnetar-strength fields in 
the heating region due to quantum effects restricting the elec- 
tron(positron) phase space (Lai & Qian 1998; Duan & Qian 
2004). We neglect the effects that high B have on q v and defer 
study of these effects to future work. In addition, for strong 
surface magnetic field strengths, heating via the dissipation 
of convectively-excited MHD waves may become important 
(Suzuki & Nagataki 2005). We assess the importance of wave 



heating and momentum deposition in ^4.3.11 

In this work we include gravitational redshifts, radial 
Doppler shifts, and modifications to the "effective solid an- 
gle" (and hence local neutrino flux) presented by the neu- 
trinosphere in the curved spacetime. The latter effect is de- 
scribed in Salmonson & Wilson (1999), while the Doppler 
and redshifts can be combined into the simple, approximate 
prescription relevant for all neutrino species: 

(0={^ D rHO\r=R v , d3) 

where 

<fc = a(R u )/a(r) , <fe = j(l-v r /c), (14) 

a(r) = yf\-(2GM/c 2 r), 7" 1 = ^l-O^ + v^/c 2 , and we 

have assumed that at radii where the equatorial flow becomes 
mildly relativistic, typical neutrinos will primarily be mov- 
ing radially. We emphasize that while we include neutrino 
gravitational redshifts in calculating heating rates we calcu- 
late the wind dynamics in Newtonian gravity. Including the 
effects of the deeper general-relativistic (GR) potential of a 
Schwarzschild metric lowers M and increases S a , Y", and the 
asymptotic wind speed V (Fuller & Qian 1996; Cardall & 
Fuller 1997). 

Our model's equation of state (EOS) includes contributions 
from photon radiation, ideal nucleons, and relativistic, degen- 
erate electrons and positrons. Non-relativistic nucleons gen- 
erally dominate the EOS near the PNS surface, but within 
a density scale height above R u the flow becomes radiation- 
dominated. As with the heating/cooling rates, high magnetic 
field effects on the EOS are ignored. 

2.3. Numerical Method 

In the steady-state Weber-Davis wind, three critical points 
occur in the radial momentum equation at radii where the out- 
flow velocity matches the local phase speed of infinitesimal 
fluid disturbances. The steady-state eigenvalues M,£, and 
B (eqs. J6), 0, and HI]) are fixed by the requirement that 
the solution pass smoothly through the slow-magnetosonic, 
Alfvenic, and fast-magnetosonic points. Physically, we 
choose these solutions over sub-magnetosonic "breezes" be- 
cause we assume the SN shock or fallback pressures at large 
radii are insufficient to stifle the strong ram-pressures of the 
wind. However, fallback at early times is not well-understood 
because it depends sensitively on the mechanism for launch- 
ing the SN shock (Chevalier 1989; Woosley & Weaver 1995) 
and thus this issue deserves further attention. Although the 
PNS radius, rotation rate, and neutrino luminosity evolve in 
time, for realistic wind conditions the timescale required for 
any of the MHD wavemodes to traverse all critical points is 
always much shorter than the timescale over which the wind 
characteristics appreciably change (e.g., tkh or the spin-down 
timescale tj). For this reason a time-series of steady-state so- 
lutions is generally sufficient to accurately model the wind 
during all phases of the PNS evolution. However, precisely 
because all physical solutions must pass through each criti- 
cal point, in the time-independent formulation of this problem 
boundary conditions must be placed on the wind solution at 
these locations. To avoid this numerically complicated singu- 
larity structure, we have instead solved the more complete, 
time-dependent version of the problem using the 6th order 
space/3rd order time, "inhomogeneous" 2N-RK3 scheme of 
Brandenburg (2001). 



5 



Our code evolves the variables (p, r, v r , v^fi^,). The value 
of the PNS mass M, neutrinosphere radius R„, magnetic flux 
= B V R 2 V , stellar rotation rate fl, and neutrino luminosity 
L v are the parameters that uniquely identify a wind solution. 
We use M = 1.4 M Q and R v = 10 km in all of our calcula- 
tions. Because our code is intrinsically non-conservative, we 
use the constancy of M, C, and I and equation ([Toll's con- 
straint on B as independent checks on the code's numerical 
accuracy. For numerical stability, an artificial viscosity of the 
form vV 2 is included in the evolution of each variable, where 
v is an appropriately-scaled kinematic viscosity (e.g., Bran- 
denburg 2001). 

We chose the location of the outer boundary, generally at 
r 1=3 1000 km, as a compromise between the run-time to reach 
steady-state and the desire to minimize the effects of artifi- 
cially forcing the fast point on the computational grid (see dis- 
cussion in i)2.3.1t . We space the radial grid logarithmically, 
choosing the number of grid points (generally 500-2000) 
to obtain the desired level of conservation of M,C,T, and B 
while simultaneously maintaining large enough artificial vis- 
cosity to maintain code stability. With sufficient resolution 
and low enough viscosity the code shows radial conservation 
of all eigenvalues to < 1% across the entire grid, although 
we did not require this level of conservation for all solutions 
so that we could efficiently explore the parameter space of 
wind properties. The mass loss rate was the most difficult 
eigenvalue to conserve yet is accurate (relative to its fully con- 
verged value) to at least ~ 10% for all solutions presented in 
this paper. 

2.3.1. Boundary Conditions 

The azimuthal speed at R u is set to enforce v$ v = 
R v £l + v r , v B$ M /B v , where a subscript v denotes evaluation at 
R v and where SI is the stellar rotation rate that defines the an- 
gular speed of the rotating frame in which the surface electric 
field vanishes (MacGregor & Pizzo 1983). In all cases we 
consider, ,„ ~ R„VL. The temperature at the PNS surface 
T{R U ) (generally ps 5 MeV) is set by requiring that the net 
heating rate q v vanish; this assumes the inner atmosphere is 
in LTE (Burrows & Mazurek 1982). Given our assumption 
that matter at the inner grid point maintains LTE, we fix the 
density at R v to be ~ 10 12 g cm" 3 so that the neutrino optical 
depth t v at the PNS surface reaches ~ 2/3, thereby defining 
a neutrinosphere. For slow rotation we find that the solution 
outside the inner few scale heights remains relatively insensi- 
tive to our choice for the inner density, although we find that 
M depends somewhat sensitively on p{R y ) as the PNS rotation 
rate increases approaching break-up. 

If all three critical points are captured on the numerical grid, 
the outer boundary conditions are not in causal contact with 
the interior wind and will have no effect on its steady-state 
eigenvalues. However, as the temperature of the wind de- 
clines, the fast magnetosonic point moves to very large radii; 
in fact, as the sound speed c s — * the fast point formally 
approaches infinity (Michel 1969). For this reason, the fast 
point is difficult to keep on the computational grid. Solutions 
without the fast point captured on the grid are sensitive to the 
outer boundary condition, with different choices altering, for 
instance, the spin-down rate. Therefore, to artificially force 
the fast point on the grid we increase the outer radial velocity 
boundary-condition until the fast point is captured. Otherwise 
equivalent solutions with the fast point naturally located and 



artificially placed on the grid were compared in several cases; 
we found that although the velocity structure changes at radii 
far outside the Alfven radius, our imposed boundary condi- 
tion had little effect on the eigenvalues of the problem and the 
correct asymptotic speed was obtained (albeit prematurely in 
radius). Since the eigenvalues uniquely determine the steady- 
state solution, this technique, when necessary, was a useful 
expedient to obtain the desired transmagnetosonic solution. 

3. RESULTS 

FigureQJsummarizes the physical regimes of PNS winds as 
a function of Lp e and rotational period P for a representative 
strongly magnetized PNS: B v = 2.5 x 10 14 G, M = 1.4 M©, 
and R v = 10 km. A cooling PNS of fixed surface field will 
traverse a path from high to low L v in this diagram, reach- 
ing Lp fi 5i ~ 0.1 at t = tkh (see eq. |28l ). If the spindown 
timescale tj is less than t kh , the PNS evolves to higher P dur- 
ing tkh, but otherwise, the PNS evolves from higher to lower 
L v at roughly constant P (see Tables [TJ and [2] for represen- 
tative rj). The regions in F igure [TJ correspond to the differ- 
ent wind phases outlined in j HUJ (1) a thermally-driven wind 
at high L v and long P; (2) a non-relativistic, magnetically- 
driven wind at high L v and short P; and (3) a relativistic, 
magnetically-driven wind at low L v and short P. In addition 
to these different regimes, Figure QJillustrates the range of ro- 
tation periods for which M is enhanced by centrifugal flinging 
(roughly P < 2 - 3 ms; see also Fig. [6J and for which Td yn , the 
dynamical time at T = 0.5 MeV (eq. l35l ). and the asymp- 
totic wind entropy S a are altered from their non-rotating, non- 
magnetized values (which has important consequences for r- 
process nucleosynthesis; see ^4. 31 . In this section, we present 
and discuss the properties of solutions for a range of parame- 
ters (B„,P,Lu e ) that span each of the non-relativistic regimes 
illustrated in Figure [TJ Some of the properties of wind solu- 
tions at Ly e = 8 x 10 51 ergs s" 1 , which corresponds to a rela- 
tively early stage in the PNS cooling evolution, are given in 
Table [TJ Table [2] compares the properties of wind solutions 
with B v = 2.5 x 10 15 G at two different neutrino luminosities 
(L Pe = 8 x 10 51 and 3.5 x 10 51 ergs s" 1 ). 

3.1. Thermally-Driven Winds 

The thermally-driven region in FigureQJcorresponds to con- 
ditions under which the PNS outflow is driven primarily by 
neutrino heating; the magnetic field and the rotation rate are 
unimportant in either accelerating or setting the mass loss rate 
of the wind. Figure [2] shows the velocity structure of such 
an effectively non-rotating, non-magnetized (NRNM) solu- 
tion for Lp e)5 i = 8, B v = 10 13 G, and tt = 50 s" 1 (P ~ 126 ms). 
Notice that the Alfven radius is relatively close to the PNS 
surface (Ra ~ 20 km) and that the sonic point (correspond- 
ing to the fast point in the NRNM limit) is at a much larger 
radius (R s w 750 km, approximately the "Parker radius", 
R p ~ GM/(v a ) 2 , of an equivalent polytropic wind, where v" is 
the asymptotic wind speed). Although the magnetic field and 
rotation rate are low enough that they have no effect on the 
wind energetics, the Alfven radius Ra is still above the PNS 
surface. Angular momentum loss is thus enhanced by a factor 
of (Ra/R„) 2 « 4 over an unmagnetized wind. In this respect, 
PNS winds such as that shown in Figure [2] are analogous to 
the solar wind (which is also primarily thermally-driven, but 
hasR A >R Q ). 
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FIG. 1 . — The regimes of PNS winds in the space of La e (= Lj} e 51 X 10 51 ergs s ') and rotation period P for a representative strongly magnetized PNS: 
B v = 2.5 X lO 14 G, R v = 10 km, and M= 1.4 Mq. A monopole field geometry is assumed. A cooling PNS will evolve from high Lc l to low La e in this diagram, 
with Lv e .51 ~ 10 at a time tQ ~ Is following the launch of the SN shock to L^ c 51 ~ 0.1 by the end of the Kelvin-Helmholtz phase (trh ~ 10— 100s after 
explosion). The solid line (eq. 1181 ) shows the boundary between PNS winds that are primarily thermally-driven by neutrino heating (high L v \ slow rotation) 
and winds that are primarily magn eto-centrifugally driven (lo w L v \ rapid rotation); for the latter, the wind power E a is enhanced (see eq. 1201 ) and the dynamical 
time Tdyn is reduced (see Fig. 1121 . The dotted line (<r = 1; eq. 1271 ) shows the boundary between non-relativistic and relativistic magnetically-driven winds. For 
sufficiently rapid rotation (P < 2 — 3 ms; t he d ot-dashed line), the mass-loss from the PNS (M) is enhanced because of centrifugal flinging (eq. 1241 ) and the 
asymptotic wind entropy S" is reduced (eq. 1371 ) because matter moves more rapidly through the region of significant neutrino heating. The dynamical time and 
entropy are important for nucleosynthesis in the wind (94.31 . For B„ > 2.5 X 10 14 G the thermally-driven region shrinks (to longer P and higher L„ e ; eq. 1181 ) 
and the relativistic region expands (to higher L^. ; eq. |27| ). For R v > 10 km, as will occur at early times when the PNS is still contracting, the wind is likely 
to be thermally driven. This figure illustrates the wide range of conditions under which PNS winds will be magnetically-driven, although it should be cautioned 

that the surface dipole field B„ v associated with the monopole field B v scales as 6^ lp oc B V P (see eq. 1301 ). which means that, for large P, the true dipole field 
appropriate to this diagram is much greater than the monopole value of B v = 2.5 X 10 14 G. 



Many of the relevant results for NRNM PNS winds (such as 
Fig.|2|i are approximated analytically and verified numerically 

in QW. QW show that M oc L% 3 (e Pe ) 10 / 3 (their eqs. [58a,b]); 
hence, if L v cx (e„) 4 , as we have assumed, QW find that 
the neutrino-driven mass loss rate is approximately given by 
Mnrnm ^3x lO- 4 (Lp l ,,5i/8) 2 - 5 M s _1 , where U e = L 9efi \ x 
10 51 ergs s" 1 . Our calculations find that Mnrnm oc L 2 ' 5 as well, 
but with a normalization lower than that of QW: 



indeed, from Figure |2]we find an asymptotic speed v NRNM w 
0.06 c at L Sti 5i = 8. T01 found that v NRNM ~ 0.1 c (Lp.^i/8) - 3 
in GR at high L v . Although we find that v NRNM scales in ap- 
proximately the same way with L p ,, our asymptotic speeds 
are lower than those obtained by T01 primarily because we 
have used a more shallow, Newtonian gravitational potential. 
In NRNM winds the asymptotic wind power is entirely gas 
kinetic energy: £ N rnm - ( 1 / 2 Wnrnm(v' N rnm) 2 ; hence, from 



1.4 x lO" 4 (L Pei51 /8) 2 ' 5 M s- 



(15) 



our results for M N rnm and v NRNM we find 



primarily because we have included neutrino redshifts in our 
heating rates. Somewhat coincidentally, T01 found a result 
similar to equation (fl~5T > from calculations incorporating GR. 

Because neutrino heating is so concentrated near the PNS 
surface, NRNM winds are barely unbound in comparison to 
the PNS escape speed (non-relativistically, VescC^) ~ 0.64 c); 



F a 

£ NRNM 



4x 10 4/ (L Pi „5i/8r 2 ergss 



(16) 



Since the time spent at L Se 51 ~ 8 is only ~ 1 s (see §4TK tne 
total energy extracted during the Kelvin-Helmholtz epoch in 
a NRNM wind is - 10 47 - 10 48 ergs. 

3.2. Magnetically-Driven Winds 
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FIG. 2. — Velocity profile for a thermally-driven wind with = 8, 

B„ = 10 13 G, and f2 = 50 s _I (P 130 ms). The variables v, , vu, Va> and Vf are 
the radial, azimuthal, Alfven, and fast magnetosonic speeds, respectively; the 
fast( Alfven) speed is also approximately the adiabatic sound(slow) speed for 
thermally-driven winds. This solution has Ms; 1.4x IO^Mq s"',ffR!3x 
10~ 8 , E" ~ 4 X 10 47 ergs s _1 , and t, ~ 880 s (see Table[T). The Michel speed 
for this solution is vm = <t'^ 3 c ~ 0.003 c, which is less than the thermally- 
driven asymptotic speed actually obtained (f» 0.06 c); hence, the magnetic 
field and rotation rate have no significant effect on the acceleration of the 
wind. 
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3. — Velocity profile for a magneto-centrifugally-driven wind with 
£i/ e ,5I = 8,B„ = 10 15 G, and O = 5000 s" 1 (P ~ 1.26 ms); the variables v r , va, 
va, Vf, and v 8 are the radial, azimuthal, Alfven, fast, and slow magnetosonic 
speeds, respectively. This solution has M ~ 3.0 X 10~ 3 Mq s~', a ~ 0.16, 
rj ~ 9 s, and Z?" ~ 2.3 X 10 51 ergs s -1 , Comparison plots of v r and 
(dashed lines) are for a 7 = 1.15 polytropic wind with similar M, Q, B v , and 
inner temperature. Notice that the radius of the slow point (approximately the 
sonic point; where v,- = v s ) is very close to t he value expected in the magneto- 
centrifugal limit: R s cf ss 19.6 km (eq. fl9l ). 

For a PNS with a given neutrino luminosity, larger rotation 
rates and magnetic field strengths lead to additional accelera- 
tion in the outer, supersonic portions of the wind. If the rota- 
tion rate and magnetic field are high enough, they will dom- 
inate the wind acceleration at large radii. This is the "Fast 
Magnetic Rotator" (FMR) limit, using the stellar-wind termi- 
nology of Belcher & MacGregor (1976); see also Lamers & 
Cassinelli (1999). An approximate criteria for this limit is 
that the magnetically-driven asymptotic speed, given roughly 



by the Michel speed v M = i&J&P/M) 1 ! 3 = a^ 3 c (Michel 
1969), exceeds the asymptotic speed obtained if the wind 
were entirely thermally-driven by neutrino heating (v^rnm)- 
Using our result for vj^ RNM , the wind will be in the FMR limit 
for magnetizations 



a>2x 10- 4 (L R „5i/8) ' 9 = a FM R 



(17) 



Using equation ( fT3T > to relate L„ r and M, our calculations im- 
ply that the PNS wind is magnetically-driven below the criti- 
cal rotation period 



/ J FM R^15(L R ,.5i/8)- I ' 7 Bi 4 ms, 



(18) 



where B v = B„x 10" G. The P = Pfmr boundary in Figure Q] 
is marked by a solid line. The magnetically-driven regime en- 
compasses a large range of PNS parameter space and hence 
generically describes most of a strongly magnetized PNS's 
evolution. In contrast, the wind from a relatively weakly mag- 
netized PNS will only be dominated by magneto-centrifugal 
forces late in the cooling epoch. The conditions necessary 
for the magnetically-driven phase to dominate the total en- 
ergy and mass loss during the Kelvin-Helmholtz phase will 
be discussed further in ^4. 1 1 



Figure [3] shows the velocity structure of a magnetically- 
driven wind from a PNS with Lp e . 5l = 8, B u = 10 15 G and 
fl = 5000 s" 1 (P ~ 1.3 ms). As the profile of in Figure 
[3] indicates, the wind corotates to w 25 km, which is far in- 
side Ra w 46 km because the magnetic field carries a sig- 
nificant fraction of the angular momentum. In addition, be- 
cause the wind is magnetically-driven, the wind speed at large 
radii is almost an order of magnitude larger than in a NRNM 
wind: v r = 0.54 c w \>m obtains at the outer grid point. The 
sonic point of the wind (corresponding to the slow point in 
the FMR limit) is now inside Ra, less than one stellar radius 
off the surface; this is expected because analytic considera- 
tions show that with increasing il, the location of the sonic 
radius R s decreases from a value of order the Parker radius to 
a value independent of the local thermodynamics (see Lamers 
& Cassinelli 1999) 5 : 



(19) 



(GM/fi 2 ) 1/3 ~m£{ 3 km. 

For comparison with our solution, Figure [3] shows the veloc- 
ity structure of an adiabatic wind (7 = 1.15) with approxi- 
mately the same M, fi, B v , and surface temperature as our 
neutrino-heated wind. The adiabatic solution agrees well with 
the neutrino-heated solution because, although M is primarily 
set by L v , once M is specified the velocity structure of the 
magnetically-driven wind becomes relatively independent of 
the details of the neutrino microphysics. This agreement im- 
plies that we can accurately map our ID neutrino-heated cal- 
culations onto multi-dimensional polytropic calculations that 
employ similar boundary conditions and a similar effective 
adiabatic index (see $4.1h . 



For FMR winds, the Michel speed obtains at large radii and 
the asymptotic wind power is therefore enhanced relative to 
equation ( fl6b : 



E a ¥MR = MB a ~(3/2)Mv 2 M : 
10 50 < 3 Ml/V^/ 3 ergss- 



(20) 



5 The subscript "cf", here and belo w, stands for "centrifugal" and relates 
to the limit described by equation )25l . 
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FIG. 4. — rj = QJ/E" is the r atio of spin-down power lost by the PNS to 
the asymptotic wind power (eq. 1121 ). This figure shows rj as a function of 
rotation period P for L& e = 8 X 10 51 ergs s and monopole magnetic field 
strengths B v = 10 13 G (cross), 10 14 G (triangle), 2.5 X 10 14 G (diamond), 
and 1()' 5 G (asterisk). For the highly magnetized solutions (with R A well 
off the surface) nearly all of the extracted rotational energy escapes to large 
radii (rj ~ 1), but for low B v and high Q most of the rotational energy is 
used to unbind the wind and hence T] 1. For slowly rotating, thermally- 
driven winds (low B„ and f2) rj <C 1 because the rotational power lost by the 
PNS is insignificant in comparison to the thermal energy supplied by neutrino 
heating. 

where M = M_ 3 x 1(T 3 MoS" 1 and B" = E a /M = (v 3 M /v a ) + 
(V) 2 /2 ~ (3/2)v| f is the Bernoulli integral at large radii in the 
FMR limit, and is composed of 2/3 magnetic and 1/3 kinetic 
energy. 

To calculate the angular momentum lost by the PNS, 
we note that for any super-Alfvenic outflow, equation (|7), 
equation ©, and conservation of magnetic flux require 
that the conserved specific angular momentum obey C = 
£IR\, where Ra is defined by the position where the radial 
outflow speed matches the radial Alfven speed: Vi-(Ra) = 
B r (RA) /V^PiRA) = va- We estimate the location of the 
Alfven point in terms of rj = HJ/E a = Vl 2 R 2 A /B a defined in 
equation ( fT2b : 

R 2 A n 2 ~(l/2) n v 2 M , (21) 

so that 6 

fl A | FMR * 1 1 v 1/2 B]{ 3 MZl /3 P^ 3 km. (22) 

Figure [4] shows rj for our wind solutions with Lp e ,5i = 8 for 
several surface magnetic field strengths (see also Table[T|i. Re- 
call that 1 /rj represents the fraction of the extracted rotational 
energy from the PNS available to energize the surrounding en- 
vironment. In the limit that vu > v esc , we find that rj « 1 and 
almost all of the rotational energy lost by the PNS emerges as 
asymptotic wind power; this limit is normally assumed when 
considering magnetized stellar spin-down. However, as the 
low B v solutions in Figure [4] illustrate, winds with short ro- 
tation periods and Vnrnm < vm < v esc can be magnetically- 
driven and yet have r\ S> 1 . The primary reason for this is that 
at high f2 the neutrino heating rate per unit mass is signifi- 
cantly reduced below its NRNM value because centrifugally- 
accelerated matter spends less time in the region where neu- 
trino heating is important. Because the wind absorbs less of 

6 TCQ took ttR A = v r (R A ) = \'m and assumed rj = 1. We find that in some 
circumstances the rj = 1 assumption is not applicable, even in the FMR limit 
(see Fig. [4}. 




P(ms) 

FIG. 5. — Spin-down timescale tj = U/Q as a function of rotation period P 
ioihu e = 8 X 10 51 ergs s~' and 4 monopole surface magnetic field strengths: 
B v = 10 13 G (cross), 10 14 G (triangle), 2.5 X 10 14 G (diamond), and 10 15 
G (asterisk). The decrease in tj for rapid rotation is due to the exponential 

enhancement in M for P < 2-3 ms (eq. 1241 ). The surface dipole field B'i' p 

ass ociat ed with the effective monopole field B v scales as bJ p oc B v P (see 
eq. 1301 ). which implies that the true dipole field appropriate to this figure 
can be much greater than B v for large P. Note that an approximate anal ytic 
expression for tj in t he magnetically-driven limit (P < Pfmr; eq- 1181 ) is 
given in equation (23). 

the neutrino energy, the magnetic field becomes more impor- 
tant for unbinding the matter from the PNS. Consequently, 
only a fraction of the rotational energy extracted at the PNS 
surface reaches large radii. 

In the limit of thermally-driven solutions with very low 
B v (even lower than in Fig. |2] such that the Alfven ra- 
dius is interior to the stellar surface), R v is the lever arm 
for angular momentum loss; thus C = flR 2 and hence r\ = 
n 2 RlMNRNM/E^ RNM sa 24(L Pf , 51 /8)-°- 6 J p-2 (using eqs. E) 
and fl6l ). For P>lms thermally-driven winds have 77 <C 1, 
which explains why 77 decreases rapidly for solutions with 
large P in Figure |4] Physically, this is because for slow ro- 
tation rates the rotational power lost by the PNS is insignifi- 
cant in comparison to the thermal energy supplied by neutrino 
heating. 

The rate at which angular momentum is extracted from the 
PNS is J = IQ = QR A M, where J = IQ is the angular momen- 
tum of the PNS and / ~ (2/5)MRl is the PNS moment of in- 
ertia. Hence, given the Alfven radius from equation d22b . the 
spin-down time of the PNS (77 = SI/ SI) in the non-relativistic, 
magnetically-driven limit is found to be 

tj | FM r =s 44077- 1 B'^MZfp^ 3 s. (23) 

This result shows explicitly that while increasing the mass- 
loading places a greater strain on the field lines (Ra oc M" 1 / 3 ) 
and hence reduces the net loss of angular momentum per gram 
(C = flR 2 oc M" 2 / 3 ), the additional mass loss carries enough 
total angular momentum to increase the overall spin-down 
rate. The high mass loss accompanying the Kelvin-Helmholtz 
epoch can therefore efficiently extract the rotational energy of 
the PNS. 

Figure [5] shows tj calculated directly from our wind solu- 
tions as a function of P at Lp c .5i = 8 for several magnetic field 
strengths (see also TableQ]). Our numerical results agree well 
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FIG. 6. — Mass loss rate M as a function of the rotation rate CI at 
Lp e = 8 X 10 51 ergs s" 1 for B v = 10 13 G (cross), fi„ = 10 14 G (triangle), 
2.5 X 10 14 G (diamond), and 10 15 G (asterisk); also shown is M(VL) for 
ho e = 3.5 X 10 51 ergs s" 1 and B v = 2.5 X 10 14 G (dotted, diamond). M in- 
creases with increasing Q and B v because centrifugal support expands the hy- 
drostatic atmosphere (see Fig. [7). For sufficiently large B v (> B c f; eq. 1261 ). 
however, M(f2) no longer increases with increasing B v because the wind 
corotates past the sonic point. An a pproximate fit to the numerical results 
in this limit is given by equation )24t . 

with the analytic estimates in equation (l23l for winds that are 
in the FMR limit (i.e., P < Prvir; eq. fTST ), The sharp decline 
in tj at short P is due to the fact that tj oc M -1 / 3 and that M 
is enhanced by centrifugal flinging for rapid rotation. Note, 
however, that only for B v > 10 15 G is tj ~ tkh ~ 10 s for 
a millisecond rotator. For B v = 10 13 G and large P the solu- 
tions are primarily thermally-driven and tj is independent of 
P; this occurs because both the Alfven radius Ra and mass- 
loss rate M (and hence rj) are independent of rotation rate for 
thermally-driven winds. 

The Alfven radius and spin-down times calculated in equa- 
tions d22i > and d23l depend on the mass loss rate M from the 
PNS, which itself depends on the PNS's rotation rate and 
magnetic field strength. Figure [6] shows our determination of 
M as a function of f2 for field strengths B v = 10 13 , 10 14 , 2.5 x 
10 14 , and 10 15 G at L p „ 5 i = 8, and B u = 2.5 x 10 14 G at 
Lv e ,s\ =3.5. The mass loss rate increases rapidly with rota- 
tion for P < 2-3 ms and also increases with B v , though it sat- 
urates for the largest magnetic field strengths, as can be seen 
by comparing the B v = 2.5 x 10 14 G and 10 15 G solutions. For 
sufficiently large B v , such that Ra > R s , we find empirically 
that M is given by 

M ~ M NRNM exp [(n/n c{ ) 2 ] ee Mrf, (24) 

where Mnrnm is the mass-loss rate for NRNM winds 
(eq. El) and O rf « 2700(L Pf , 5 i/8) ' 08 s" 1 . 

The enhanced mass loss shown in Figure [6] is due to the 
effect of strong magnetic fields and rapid rotation on the sub- 
sonic, hydrostatic structure of PNS winds. Figures [7] and [8] 
show the density and temperature profiles, respectively, for 
winds with B v = 10 15 G and Lp ei 5i = 8 at several rotation rates. 
For the most rapidly rotating solution (solid line), Figure [7] 
shows that centrifugal support is sufficient to expand the scale 
height of the hydrostatic atmosphere at small radii, resulting 
in the much higher mass loss rates seen in Figure [6] Analyti- 
cally, we expect the centrifugal support to be important when 
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FIG. 7. — The density profiles in the inner ~ 100 km for magnetically- 
driven wind solutions at L^. 5 i = 8, B v = 10 15 G, Q. = 7000 s _1 (solid), 5000 
s" 1 (dotted), 2000 s" 1 (dashed), and 500 s~' (dot-dashed). Tr iangles, crosses, 
and diamonds mark the slow point, the radius where T = 0.5 MeV, and the 
Alfen radius, respectively. For all of the solutions, the supersonic portion of 
the density profile is altered from that of a thermally-driven solution due to 
magnetic acceleration. For solutions with £1 > 3000 s~' the density scale 
height in the subsonic portion of the wind (interior to the triangle) is signifi- 
cantly larger due to centrifugal support. This enhances the mass loss rate as 
shown in Figure[6] 

Rl,£1 > c S) „, where c s „ is the sound speed at the PNS neutri- 
nosphere. We find that the inner sound speed depends only 
weakly on the neutrino luminosity: c s „ ps 0.12(Lp e) 5i/8) 
c, and thus that mass loss is enhanced for 

fi > 3600(L Pf ,5i /8) 08 s _1 . (25) 

This region is denoted "M Enhanced" in Figure [T] Equation 
( |25] | is in good agreement with our numerically determined 
value of fl c f defined in equation d24l . 

The enhancement of M implied by equations (f24b and d25l l 
does not explicitly depend on the magnetic field strength. 
However, centrifugal support of the PNS atmosphere only oc- 
curs if the field can sustain corotation out to the sonic radius 
R s ; otherwise rotation has much less of an effect on the mass 
loss rate (see the solutions with low B v in Fig. [6j. The re- 
quirement of corotation out to R s can be written as R A > R s C f 
(eq. fl9l ), which in turn implies B v > B ct , where 

B ci ~ 2 x 10 14 ^Wfp^ 2 G. (26) 

For B u > B c f , the mass loss from the PNS no longer increases 
with increasing B v (see Fig. |6]l because the wind already 
co-rotates out to the sonic radius where M is set. In deriv- 
ing B c f in equation d26| i, we have used equation ( l22l for Ra 
because, under most conditions, a wind that is centrifugally- 
supported will automatically be in the FMR limit, although 
the converse is not necessarily true. We note that even rela- 
tively weakly-magnetized PNS's that are rapidly rotating will 
experience some degree of enhanced mass loss. For instance, 
for 17 = 7000 s , even the lowest field strength solution in 
Figure[6](B„ = 10 13 G) has a mass-loss rate almost an order of 
magnitude larger than its NRNM value. 

The non-relativistic calculations we have presented here are 
only applicable for magnetizations a < 1 ; this is equivalent to 
requiring Ra </?l = c/Vt, the radius of the light cylinder. For 
a > 1 the PNS wind becomes relativistic and its spin-down 
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FIG. 8. — The temperature profiles in the inner ~ 100 km for magnetically- 
driven wind solutions at L^, 5 i = 8, B v = 10 15 G, fi = 7000 s _1 (solid), 5000 
s" 1 (dotted), 2000 s" 1 (dashed), and 500 s" 1 (dot-dashed). Tr iangles, crosses, 
and diamonds mark the slow point, the radius where T = 0.5 MeV, and the 
Alfen radius, respectively. Notice that in all but the most rapidly-rotating case 
the T = 0.5 MeV radius is located between the slow point and the Alfven ra- 
dius, which implies that there is significant mag neti c acceleration of the wind 
at the T = 0.5 MeV radius. As shown in Figure [T2l this significantly reduces 
the dynamical time at T = 0.5 MeV, making magnetically-driven PNS winds 
more favo rable for r-process nucleosynthesis than thermally-driven winds 
(see S3). 

properties will change (see ^4.1 I f or a discussion). Using equa- 
tion d24b for M (i.e., assuming B v > B C f), the magnetization is 
given by 

a ~ 0.05^ 4J p-2(L Pe , 51 /8)- 2 - 5 e X p[-5.4 J p-2(L Pe . 51 /8)- 016 ]. 

(27) 

The (7=1 boundary is denoted by a dotted line in Figure 
Q] PNSs with a > 1 and Lp,,5i > 0.1 will experience a rel- 
ativistic phase accompanied by significant mass loss; this 
mass loss keeps the wind mildly relativistic, in contrast to the 
much higher a spin-down that will commence following tkh 
(L Pei5 i «0.1). 

4. APPLICATIONS AND DISCUSSION 

4.1. Magnetized PNS Evolution 

With our numerical results in hand that sample a wide range 
of PNS wind conditions, we can begin to address the time 
evolution of a cooling, magnetized PNS. In the early stages 
following the launch of SN shock the PNS is likely hot and 
inflated, with a radius exceeding the value of = 10 km that 
we have assumed in all of the calculations presented in this 
paper. This early phase is likely to be thermally-driven for 
all but the most highly-magnetized proto-magnetars, and, if a 
dynamo is at work, the large-scale field itself may still be am- 
plifying during this phase (Thompson & Duncan 1993). Us- 
ing the collapse calculations of Bruenn, De Nisco, & Mezza- 
cappa (2001) (from a 15 M Q progenitor of Woosley & Weaver 
1995) T01 fit an approximate functional form to the PNS ra- 
dial contraction: R v oc f -1 / 3 such that R v (l s) ~ 15 km and 
R v (2 s) ~ 12 km. The recent SN simulations of Buras et 
al. (2003, 2006) with Boltzmann neutrino transport find a sim- 
ilar neutrinosphere radius at t ~ Is after bounce. 

We do not attempt to address the uncertainties in early-time 
PNS cooling calculations, especially in the presence of large 



fields and rapid rotation. Rather, we assume that the PNS 
has cooled to its final radius and completely established its 
global field by a time to ~ 1 s following core-collapse; we can 
then use the calculations presented in this paper to investigate 
the subsequent evolution of the PNS. Following ?o we assume 
a simplified PNS cooling evolution (similar to that used in 
TCQ, motivated by Figure 14 of Pons et al. 1999): 

Lp, j5 l(0=-tfl( — ) ■■t <t<TKH 

\tkhJ 

Lu e ,5i(t) = Lo exp[-(f-T KH )/TkH] : * > t kh , (28) 

where, for definitiveness in what follows, we take Lo = 0.2, 
to = 1 s, tkh = 40 s, and 6 = 1 . This cooling evolution is ap- 
proximate because magnetar-strength fields and rapid rotation 
could alter r KH or the form of the cooling profile (e.g., 5) by 
affecting the neutrino opacity or the dynamics of the contrac- 
tion itself (e.g., Villain et al. 2004; Dessart et al. 2006). For 
instance, in ID collapse calculations with rotation, Thomp- 
son, Quataert, & Burrows (2005) found that for P ~ 1 ms, the 
total neutrino luminosity at t ~ 0.6 s after bounce is ~ 50% 
smaller than in a non-rotating PNS. 

The dominant uncertainty in applying our results to magne- 
tized PNS evolution is that we have assumed a monopole field 
geometry. To relate our results to more realistic dipole sim- 
ulations, we use the recent axisymmetric, relativistic MHD 
simulations of B06, who simulate neutron star spin-down for 
a « 0.3-20. B06 show that the energy and angular momen- 
tum loss rates from aligned dipole spin-down can be described 
accurately by monopole formulae provided they are normal- 
ized to just the open magnetic flux; for instance, we can ac- 
curately apply our results for tj (eq. 11231 ) with a suitable 
renormalization of B v . 

To apply the results of B06 we need to estimate the open 
magnetic flux in PNS winds. In force-free spin-down calcu- 
lations motivated by pulsars it is generally assumed that the 
radius of the last closed magnetic field line (the "Y point" 
Ry) is coincident with the light cylinder (Contopoulos et 
al. 1999; Gruzinov 2005) so that the ratio between the frac- 
tion of open magnetic flux in the dipole and monopole cases is 
Ru/Ry = Ru/Rl- This assumption is supported by force-free 
simulations (Spitkovsky 2006; McKinney 2006), which show 
that when mass-loading is completely negligible (a — > oo), 
Ry — Rl- However, B06 show that rapidly rotating, mildly 
mass-loaded MHD winds have a larger percentage of open 
magnetic flux than vacuum or force-free spin-down would im- 
ply (i.e., Ry < Ri,). From Table 4 of B06 we fit the approxi- 
mate power law 

fl y /fl L ~0.31cr - 15 (29) 

for the range a £ [0.298,17.5] and for a fixed rotation pe- 
riod of order one millisecond. Although there are uncertain- 
ties in quantitatively extrapolating B06's results, reaching the 
pure force-free limit with Ry ~ R^ appears to require a ^> 1 . 
We therefore conclude that magnetized, rapidly rotating PNS 
winds (with a G {10~ 3 , 10 3 } for t < tkh under most circum- 
stances) will typically possess excess open magnetic flux. 

Because the results of B06 for Ry / R^ cover only a relatively 
narrow portion of PNS parameter-space we must proceed with 
caution in generalizing their results to our calculations; on the 
other hand, their basic conclusion shows a weak dependence 
on a and fi, and has a solid theoretical explanation (Mestel & 
Spruit 1987). Hence, we have attempted to apply the results of 



11 



B06 to gain insight into the multi-dimensional generalization 
of our calculations, but we check at every step that our cal- 
culations are not overly sensitive to extrapolations of B06's 
results. Combining equation ( f29b with the empirical formula 
B06 provide (their eq. [25]) relating the effective monopole 
surface field B™ on to the true equatorial dipole field p we 
find 

B7*~0.6B^^ ~0.4B d ;Pp,>- (U5 , (30) 

which remains approximately valid for a substantial range in 
er, provided we enforce Ry = for a > 10 3 and keep Ry > 
Ru- 

Using equation (l30l > we substitute B™ on for B v in the re- 
sults of $3] and integrate from to to tkh using our fiducial 
cooling evolution (eq. l28l ) to obtain the total energy and 
mass loss during the Kelvin-Helmholtz epoch (E" ot and AM M , 
respectively) as a function of the fixed dipole surface field 
B^ p and for initial rotation periods Po e {1,10} ms. Al- 
though a neutron star will eventually impart any remaining ro- 
tational energy after tkh to its surroundings through an ultra- 
relativistic, pulsar-like outflow, we concentrate on the wind 
evolution prior to the end of the Kelvin-Helmholtz phase be- 
cause our calculations are primarily suited to studying mass- 
loaded spin-down and because we are interested in energy 
that can be extracted sufficiently early to affect the rapidly 
outward-propagating SN shock. In performing our calcula- 
tions when the wind is relativistic (a > 1), we continue using 
equation d24l) for M and employ the relativistic spin-down for- 
mula given by B06 (their eq. [26]): 

Z^1.5xl0 4 K?«(^) 'ergss- 1 , (31) 

evaluated for Ry/Rl given by equation (|29l >. 

Although it is possible that equation d29l ) may not be ac- 
curate far outside the parameter regime B06 considered, we 
found that re-running the calculations at fixed Ry/Rl= 1/3 
changed E° ot by, at most, a factor of a few. A similar uncer- 
tainty in our calculations is that we have used our equatorial 
M over all 4ir steradian; thus, we have probably overestimated 
AM tot by a factor of ~ 2 due to the effects of closed magnetic 
flux and because centrifugal flinging concentrates M at low 
latitudes. Although a direct comparison of the mass loss rate 
between our solutions and the dipole simulations of B06 is 
difficult, we find that the dependence of M on is similar 
between our solutions when the surface temperature in B06's 
simulations is scaled to an appropriate neutrino luminosity. 

Our estimates for E" ot and AM to t during the Kelvin- 
Helmholtz phase are presented in Figures [9] and [10] for val- 
ues of B & v V ranging from 10 12 - 10 16 G. The single domi- 
nant phase contributing the majority of the energy extracted 
for a given initial rotation rate p> is denoted by the line 
style (thermally-driven [NRNM] = dotted; non-relativistic, 
magnetically-driven [FMR] = solid; relativistic, magnetically- 
driven [REL] = 3 dot-dash). While all PNSs pass through each 
wind phase sometime during tkh (see Fig. IT}, PNSs can still 
be usefully classified into 3 types based on which wind phase 
dominates the total energy loss during t^h' 

(1) Thermally-DrivenWinds(Bii, p <10- 2 P^ s ): For low 
magnetic field strengths, Figures l9l and [TOl show that the total 



mass and energy loss are effectively at the NRNM values for a 
large range of initial periods: AM£ t RNM ~ lO- 4 (f()/ls)- L5 M , 
£tot RNM — 2 x 10 47 (fo/ls)" 2 - 2 ergs. Energy and mass loss from 
this class of PNSs is generally modest and is dominated by 
early times. Analysis of the Parkes multibeam survey sug- 
gests that half of all pulsars are born with B^ < 3 x 10 12 G 
(Vranesevic et al. 2004), implying that NRNM winds dom- 
inate the majority of neutron star births, independent of the 
birth-period distribution. Spin-down during tkh is negligi- 
ble for PNSs of this type and the supernova remnants asso- 
ciated with the production of NRNM PNSs will not be sig- 
nificantly modified by the small energy injected during the 
cooling phase. 

(2) Non - Relativistic, Magnetically - Driven Winds 
(lO^PL < B* p < 2PS ls exp[2P^]) : For ~ 10° - 10 15 G sur- 
face field strengths, Figure [9] shows that E" ot is dominated by 
a non-relativistic, magnetically-driven outflow during tkh for 
most periods between 1 and 10 ms. Note that most observed 
Galactic magnetars have field strengths in this range (Kou- 
veliotou et al. 1998). Figure [9] shows that for B — 10 15 G 
and P < 2 ms, more than 10 51 ergs can be lost to a non- 
relativistic, magnetically-driven outflow during t%h, and that, 
over a broad range of initial spin period, the energy extracted 
is many times larger than from a slowly rotating PNS. Be- 
cause non-relativistic outflows are efficiently collimated along 
the rotational axis by magnetic stresses (B06), the energy per 
unit solid angle at the pole may exceed that of the SNe, po- 
tentially altering the SNe shock's morphology and nucleosyn- 
thetic yield. 

If, for the purposes of an analytic estimate, we assume that 
Ry = Pl/3 then the total energy extracted for this class of PNS 
can be approximated as 7 

£™ R ^ 10 50 (< P ) 4/3 P-r 
x[r f 1/6 -f,; /6 ]exp[1.3P- 2 ]ergs, (32) 

where Tf = minjrKH, trel}, to and Tf are in seconds, and trel 
~ 7 s (Bij P )-°' 8 ^ s 6 exp[2.2P- 2 ] is the time after which the 
wind becomes relativistic. The weak dependence of equation 
( 1321 on tkh and to shows that the total energy extracted in the 
magnetically-dominated phase is relatively insensitive to our 
choice for the PNS thermal evolution because the energy loss 
is distributed almost equally per decade in time. For the same 
reason, P™ R is relatively insensitive to the precise time dur- 
ing the Kelvin-Helmholtz epoch at which the PNS cools to its 
final radius. 

Figure [10] shows that the total mass loss from PNSs is en- 
hanced for P < 3 ms and Bu p > 10 13 G; the total mass loss 
increases with fij p , saturating for Z?J P > 3 x 10 14 G (i.e., 
B v B c f, see eq. 11261 ). Assuming no evolution of f2, the 
integrated mass loss is approximately 

AM tot ~ 10-V 5 exp[5.4p- 2 ]M (33) 

for large B v . For Po ~ 1 ms and Po ~ 3 ms, the total mass 
extracted is thus ~ 10 2 and ~ 2 times greater, respectively, 

7 Equation )32l also assumes that <5 = 1 , that strict corotation can be main- 
tained by the magnetic field (this criteria is well-satisfied because late times 
dominate the energy release, and winds at late times have lower M and are 
easier to support magneto-centrifugally), and that Q does not evolve signif- 
icantly (i.e., tj > tkh); this is well-satisfied for field-strengths at which the 
outflow energy is indeed extracted via a non-relativistic outflow. 
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FIG. 9. — E" ot as a function of the PNS initial period Pq for fixed dipole surface fields that range from those typical of rotation-powered pulsars to those capable 
of producing hyper-energetic SNe. E" ot is the total energy carried by the PNS wind at infinity, calculated by evolving the PNS from time to = 1 s to tkh = 40 s, 

assuming the PNS cooling evolution of equation (28) and using results from the aligned dipole simulations of B06 to relate bJ p to our monopole calculations 
(see eq. |30| ). The line style denotes the wind phase that dominates the total energy loss (dotted = thermally-driven; solid = non-relativistic, magnetically-driven; 
3 dot-dash = relativistic, magnetically-driven). The dashed line at the top shows the total initial rotational energy of the PNS. For B A J V = 10 12 G neutrino-heated, 



thermally-driven outflow dominates for Pq > 1 ms, while for B^ p - 
wind during the Kelvin-Helmholtz phase. 

than from a non-rotating, non-magnetic PNS. Equation ( f33l > 
is shown with a dashed line in Figure [TOl 

(3) Relativistic, Magnetically -Driven Winds 

(Bf£ > 2P^ s exp[2P- 2 s ]) : For = 10 16 G and P < 6 ms, 
the energy loss from PNSs is dominated by a relativis- 
tic, magnetically-driven outflow during the Kelvin-Helmholtz 
phase. As Figure [9] illustrates, PNSs of this type lose ss 
10- 100% of their total rotational energy during t kh . For 
sufficiently rapid rotation such PNSs are therefore candidates 
for the central engine of hyper-energetic SNe and LGRBs (see 
^4. 2b . PNS winds of this type differ from the ultra-relativistic 
(a ^ 10 3 ), pulsar-like phase that begins once neutrino heat- 
ing completely subsides in two important ways: (1) neutrino 
heating causes significant mass loss that maintains modest a 
(< 10 2 - 10 3 ); and (2) mass-loss leads to more open magnetic 
flux (eq. 11291 ) than in the pure force-free case and thus the 
PNS spins down more rapidly (see eq. OTl ). 

4.2. Hyper-Energetic SNe and Long Duration Gamma-Ray 

Bursts 

One observational manifestation of early energy loss from 
rapidly rotating proto-magnetars may be "hyper-energetic" 



10 16 G almost the entire rotational energy of the PNS is extracted by a relativistic, magnetized 



SNe, which we define as having energy greater than a SN's 
usual ~ 10 51 ergs. SNe significantly energized by proto- 
magnetar winds are naturally asymmetric given the preferred 
direction associated with the rotation axis of the PNS, and if 
energized sufficiently early, their nucleosynthetic yield may 
be appreciably modified; additionally, in some cases proto- 
magnetar winds may provide conditions favorable for long- 
duration GRBs within just seconds of the progenitor core col- 
lapse (Thompson 1994; Wheeler et al. 2000; TCQ). In order 
to possess > 10 51 ergs of rotational energy a PNS must be 
born with a period Pq ^ 4 ms, and, as Figure [9] illustrates, 
such a large energy is extracted during the Kelvin-Helmholtz 
epoch only, for B^ p > 10 15 G, if the energetically-dominant 
form of outflow is at least mildly relativistic. Using equation 
OTb . the spin-down timescale for relativistic outflow from a 
proto-magnetar is 

r J REL ^ 350(<V^ s (^) 2 s. 

When Ry ~ Rl (f ^> 1), equation (l34l reduces to the canoni- 
cal force-free ("va cuum dipole") spin-down timescale, but as 
was discussed in fl4.ll neutrino-heated mass flux (which is 
significantly enhanced for P ~ 1 ms) maintains modest a at 



(34) 
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FIG. 10. — The total mass extracted via a PNS wind (AM tot ) from time 
to = 1 s to tkh = 40 s as a function of the PNS initial period Pq for dipole 
surface fields B^ p = 10 12 , 3 x 10 12 ,10 l3 ,3x 10 13 ,10 14 ,3x 10 14 , and 10 15 G 
(from bottom to top). These calculations assume the PNS cooling evolution 

of equation {28) and use the aligned dipole simulations of B06 to relate B„ p 
to our monopole calculations (see eq. |30| ). The line style denotes the wind 
phase that dominates the total energy loss (dotted = thermally-driven; solid = 
non-relativistic, magnetically-driven; 3 dot-dash = relativistic, magnetically- 
driven). The results of this figure demonstrate that for Pq < 2—3 ms and 

B* p > 10 13 G, centrifugal flinging enhances the total mass extracted from a 
PNS wind during tkh ■ An analytic approximation to AMtot for high fi„ is 
given in equation (33} and is shown with a thick dashed line. For P < 3 ms, 

AMtot is essentially the same for the B* P = 3 X 10 14 G and 10 15 G models 
because of the saturation of M at fixed Q for B > B c f . 

early times, and therefore the PNS may spin down up to an 
order of magnitude faster during the Kelvin-Helmholtz epoch 
(B06). 

For a surface dipole field typical of observed Galactic mag- 
netars (~ 10 14 - 10 15 G) the SN shock is only energized with 
> 10 51 ergs during tkh for initial periods < 1 -2 ms; how- 
ever, such a large total rotational energy (> 10 52 ergs) proba- 
bly cannot be typical of magnetar birth because even if it is not 
extracted during r KH this rotational energy will eventually be 
transferred to the surrounding environment, and observations 
of Galactic magnetar SN remnants are inconsistent with such 
a large energy (e.g., Vink & Kuiper 2006) 8 . In addition, the 
rate of hyper-energetic SNe is probably much smaller than the 
rate of magnetar births: Podsiadlowski et al. (2004) estimate 
that the hyper-energetic SNe rate is only ~ 0.01 -0.1% of the 
radio pulsar birthrate, while Woods & Thompson (2004) esti- 
mate that the Galactic magnetar birthrate is at least ~ 10% of 
the radio pulsar birthrate. 

Consider, however, a rarer class of proto-magnetar with 
rapid initial rotation (Pq ~ 1 ms) and a somewhat stronger 

global magnetic field (B d j v w 3 x 10 15 - 10 16 G; if a rapid birth 
period is indeed the cause of such a strong field, these as- 
sumptions are not independent). Figure [9] shows that winds 
from PNSs with these characteristics are dominated energeti- 
cally by at least mildly relativistic outflow. Evolving a proto- 
magnetar of this type with Pq = 1 ms using the calculations de- 
scribed in i j4.1l we find that the total energy extracted during 
the Kelvin-Helmholtz phase is 1.1 x 10 52 (1.8 x 10 52 ) ergs for 

8 One caveat to this argument is that a significant portion of the rotational 
energy could escape as gravitational waves in a time <C tj (see Stella et 
al. 2005). 



a surface dipole field strength of 3 x 10 15 (10 16 ) G; this repre- 
sents 50%(80%) of the total rotational energy of the PNS. We 
find that almost all of this energy is extracted with a < 10 3 . 
More specifically, for fi^ p = 3 x 10 I5 (10 16 ) G, we find that 
w 7 x 10 5l (l.l x 10 52 ) ergs of rotational energy is extracted 
with a < 10 in the first 5(2) s following the launch of the 
SN shock; the assumption of excess open magnetic flux over 
the pure force-free case (eq. |29l ) is therefore especially well- 
justified because wind solutions with precisely these parame- 
ters (B dip - 10 15 - 10 16 G, P ~ 1 ms, and a ~ 0.1- 10) have 
been calculated by B06. The rapid spin-down (eq. 11341 ) and 
efficient energy extraction that occurs immediately following 
the birth of a magnetar of this type may energize the SN shock 
sufficiently rapidly to enhance its 56 Ni nucleosynthetic yield 
(e.g., Nakamura et al. 2001), one of the observational signa- 
tures of hyper-energetic SNe (however, see Soderberg 2006). 
We note, however, that the ability of proto-magnetar winds to 
affect the SN nucleosynthesis is sensitive to the evolution of 
the magnetic field and radius of the PNS, since the time for 
the PNS to contract to its final radius is similar to the time 
over which energy must be extracted to affect the 56 Ni yield. 

Although a significant portion of the PNS rotational en- 
ergy emerges with a < 10, which causes enhanced spin- 
down at early times, we find that the total energy extracted 
is distributed almost uniformly in log(er) in the range a £ 
{0.1, 1000}. Indeed, for Bf v = 3 x 10 I5 (10 16 ) G we find that 

2 x 10 51 (4 x 10 51 ) ergs is extracted with 10 < a < 100 in the 
first 18(7) s, and that an additional 10 51 (3 x 10 51 ) ergs is ex- 
tracted with 100 < a < 1000 by 39(32) s following the launch 
of the SN shock. If the Poynting-Flux of this outflow can 
be efficiently converted to kinetic energy (e.g., Drenkhahn & 
Spruit 2002), then the high energy to mass density implied 
by the wind's large a will result in acceleration to a com- 
parably large asymptotic Lorentz factor. The fact that a sig- 
nificant portion of the PNS rotational energy emerges with 
a ~ 10- 1000 on a timescale rf^ - r KH ~ 10- 100 s thus 
makes the birth of proto-magnetars with Pq ~lms and B v p ~ 

3 x 10 15 - 10 16 G a viable candidate for the central engine of 
LGRBs. Indeed, it is important to note that for initial PNS 
periods that give the right energetics for LGRBs (P w 1-3 
ms) and surface magnetic fields that give the right timescale 
for LGRBs {B^ v « 3 x 10 15 - 10 16 G) the magnetization of the 
resulting proto-magnetar wind - which is not a free parame- 
ter of the problem - is consistent with the Lorentz factors in- 
ferred from LGRBs (r ~ 100; Lithwick & Sari 2001; Granot 
& Kumar 2005). One shortcoming of our calculations, how- 
ever, is that we cannot address whether proto-magnetar winds 
will have large-scale collimation, as suggested by observa- 
tions of some LGRB afterglows (e.g., "jet breaks"; Rhoads 
1997, 1999; Frail et al. 2001). As noted in the introduction, 
the collimation of the relativistic wind may depend on its in- 
teraction with less relativistic material (which is ejected ear- 
lier and is more likely to be collimated about the pole; B06) 
or the stellar mantle (e.g., Wheeler et al. 2000; Uzdensky & 
MacFadyen 2006). 

4.3. r-process Nucleosynthesis 

PNS winds are a plausible candidate for the astrophysi- 
cal source of heavy r-process nuclides (Woosley & Hoffman 
1992; Meyer et al. 1992). The chief requirement for the r- 
process in PNS winds successfully reaching the critical third 
abundance peak at A ss 195 is that the ratio of neutrons to seed 
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FIG. 1 1 . — (S") 3 /r ( iy n as a function of the rotation period P for monopole 
magnetic field strengths B v = 10 13 G (cross), 10 14 G (triangle), 2.5 X 10 14 
G (diamond), and 10 15 G (asterisk) at = 8 X 10 51 ergs s -1 , where 5° is 
asymptotic wind entropy and Tdyn is the dynamical time evaluated at T = 0.5 
MeV (eq. (35]). The ratio (S fl ) 3 /T dyn is also shown for B v = 2.5 X 10 14 G at 
La t = 3.5 X 10 51 ergs s _I (diamond, dotted). Shown with dashed lines are the 
approximate thresholds above which r-process can proceed to the third abun- 
dance peak (A fa 195), taken from the numerical study of Hoffman, Woosley, 
& Qian 1997 (their Table 5), for several 7°; notice that the threshold ac- 
tually decreases with Y* for Y" > 0.46. The dot-dashed line is (S") 3 /^™ 
calculated from the analytic expressions given by Qian & Woosley (1996) for 
NRNM winds at t,=8x 10 51 ergs s" 1 (their eqs. [48a] and [61]). This fig- 
ure highlights that strongly magnetized, rapidly rotating PNS winds produce 
conditions significantly more favorable for successful third peak r-process; 
the optimal conditions obtain for B v > 10 14 G and P ~ 2- 10 ms because 
T dyn i s reduced to a fraction of the rotation period by magneto-centrifugal 
acceleration (see Fig. 1 121 . The decrease in (5 D ) 3 /r c | yn for P < 3 ms arises 
because, for sufficiently rapid rotation, magneto-centrifugal acceleration re- 
duces the advection time of wind materi al th rough the heating region, thus 
decreasing the asymptotic entropy S°(eq. 1371 ). 

nuclei (the "neutron-to-seed" ratio) remain high until the out- 
flow cools to the point at which r-process can commence (T w 
0.1 MeV). As many previous investigations have emphasized 
(e.g., Hoffman, Woosley, & Qian 1997; hereafter HWQ), the 
primary wind properties necessary to achieve and maintain a 
large neutron-to-seed ratio are: (a) a high asymptotic wind 
entropy 9 S a , (b) a small asymptotic electron fraction Y" (large 
neutron fraction), and (c) a short dynamical timescale Td yn , 
where we follow HWQ in defining Td yn as 10 
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where the second equality only holds as an order-of- 
magnitude estimate. The dynamical time is defined at T = 0.5 
MeV because this is the radius at which a-particles, the build- 
ing blocks of seed nuclei, first form. 

HWQ present numerical calculations delineating the re- 
gions of (S a ,Td yn ,Y") parameter space necessary for nucle- 

9 More precisely the r-process requires a high wind entropy at T = 0.5 
MeV (So 5Mcv); however, entropy from neutrino heating saturates insid e the 
radius where T = 0.5 MeV so that So.sMeV — S" (in deed, a s discussed in i|2.2l 
we artificially set q v = for T < 0.5 MeV). In ! j4.3.1l we will distinguish 
between So. 5 MeV and S" when we consider the more radially-extended effects 
that wave heating can have on the wind entropy. 

10 We caution that r^ yn is sometimes defined in terms of the density scale 
height and this dynamical timescale, under the radiation-dominated and ap- 
proximately constant entropy condit ions at T ~ 0.5 MeV, is a factor of 3 
shorter than that defined in equation 135) . 



FIG. 12. — The dynamical time r^ yn at the radius where T = 0.5 MeV 
(eq. |35| ) as a function of the PNS rotation period (P) at Lp e = 8 X 10 51 ergs 
s _1 for monopole field strengths B v = 10 13 G (cross), 10 14 G (triangle), 2.5 X 
10 14 G (diamond), and 10 15 G (asterisk). The dynamical time is also shown 
for By =2.5x 10 14 GatLj7 t =3.5x 10 51 ergs s" 1 (diamond, dotted). For slow 
rotation and weak magnetic fields, Tdyn approaches a constant value (rj 20 
ms at Ly e = 8 X 10 51 ergs s _1 ). For larger B v , Td yn decreases to a fraction of 
the PNS rotation period due to magneto-centrifugal acceleration around the 
T = 0.5 MeV radius; the field required to minimize Tdyn is approximately B c f 
(eq. |26| ). For very rapid rotation (P < 2-3 ms), centrifugally-enhanced mass 
loss reduces the Alfven radius, which in turn decreases the effectiveness of 
magneto-centrifugal acceleration and increases Td™. The large reduction in 
Tdyn in strongly magnetized, rapidly rotating PNS winds provide s con ditions 
more favorable for successful r-process nucleosynthesis (see Fig. lilt . 



osynthesis to reach the third abundance peak, assuming an 
adiabatic cooling model. In general, they find that the condi- 
tion for successful r-process takes the approximate functional 
form 

S a > (r dy „) 1/3 ^, (36) 

where T is a function of only Y" (see HWQ eqs. [20a,b] 
for analytic approximations to these results). Hence, at fixed 
neutron abundance (i.e., fixed Y") a successful r-process is fa- 
vored for large values of the ratio (S ) 3 /Td yn . However, for 
a wide range of reasonable PNS properties, detailed NRNM 
wind calculations show that this ratio falls short of that re- 
quired to reach the third peak (QW; Otsuki et al. 2000; Wanajo 
et al. 2001; T01). Given that magnetar birth is relatively com- 
mon and that rapid rotation may be its key ingredient, we 
have quantified the effects that magnetar-strength fields and 
rapid rotation have on the PNS wind properties that determine 
whether third peak r-process is successful. 

Figure [TT1 shows (5 a ) 3 /r d y n from our magnetized wind so- 
lutions at L 9et5i = 8 and B u = 10°, 10 14 , 2.5 x 10 14 , and 10 15 
G as a function of the PNS rotation period; also shown is 
(S fl ) 3 /Y dyn for B v = 2.5 x 10 14 G and L Pf , 5 i = 3.5 (dashed line). 
To put our results in context, we plot the entropy threshold 
given by HWQ (their Table 5) at the lowest Td yn they consider 
(rj 5 ms; choosing such a low Td yn will be justified below) 
with a dashed line for several Y". We plot several Y" be- 
cause the entropy required for successful r-process depends 
sensitively on Y°, but a modest change to Y" in our wind cal- 
culations would not significantly alter S" or r dyn from those 
obtained with Y° = 0.4. Figure QT] shows that the presence of 
a magnetar-strength field and mildly rapid rotation moves the 
critical wind parameters almost an order-of-magnitude closer 
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to successful third peak r-process in the space of (5 fl ) 3 /Td yn - 11 
We examine the reasons for this result below. 

In rapidly rotating, strongly-magnetized PNS winds, 
centrifugal-flinging pushes matter quickly through the heat- 
ing region, which reduces S a . Our calculations find that, for 
B v > B c f (eq. |26l ), S a is approximately given by 

S a ~S a NRNM exp[-n/n s ], (37) 

where n s ~ 3500(L Pi „ 51 /8) 015 s" 1 and ^nrnm - 
70(L Po 5i/8)" ' 2 is the asymptotic entropy for NRNM 
winds in units of fc^/baryon. 

At face value, the exponential decrease in entropy implied 
by equation ( f37l > appears to stifle the chances for successful 
third peak r-process in rapidly rotating, strongly-magnetized 
PNS winds. However, because such PNS winds are mag- 
netically driven, Td yn will also decrease with increasing fi; 
hence, success for the r-process depends on the competi- 
tion between changes in Td yn and S a . As Figure [8] illustrates, 
in magnetically-driven PNS winds the T = 0.5 MeV radius 
(^o.5Mev) is generally outside the sonic point and the heating 
region (interior to which S a is set); therefore, while periods < 
2-3 ms are required to significantly affect the dynamics in the 
subsonic heating region and alter M or S a (the "M Enhanced" 
region of Fig. [TJ, Td yn is significantly reduced for more modest 
rotation rates (the much larger "Magnetically-Driven" region 
in Fig.[T). 

Figure [12] shows Td yn from our magnetized wind solutions 
at L^, 5 i = 8 and B v = 10 13 , 10 14 ,2.5 x 10 14 , and 10 15 G as 
a function of the PNS rotation period; also shown is Td yn for 
B v = 2.5 x 10 14 G and L^ S \ = 3.5 (dashed line). For rota- 
tion periods ~ 1 - 10 ms and B v > 10 15 G our solutions ob- 
tain Td yn s» 0.3-3 ms, just a fraction of the rotation period, 
while for slow rotation Td yn approaches a value w 20 ms for 
£p e ,5i = 8, similar to that derived by QW (their eq. [61]). As 
a comparison of the B v = 2.5 x 10 14 and 10 15 G solutions 
in Figure Q~2] shows, the effects of magnetic fields on Td yn 
saturate for sufficiently large fields. Indeed, the monopole 
surface field required to minimize Td yn is approximately B c f 
(eq. l26l ). Ignoring centrifugal enhancement of M and using 
equations Sl5[ and (l30l , this corresponds to a surface dipole 
field - 10 I5 (P/4ms) 3 / 2 (L Pe ,5i/8) L25 G; hence, surface fields 
typical of observed Galactic magnetars are probably sufficient 
to minimize Td yn at early times in a PNS's thermal evolution. 
Thus, if observed Galactic magnetars were born with periods 
in the range 2 ms < P < 10 ms we conclude that there is an en- 
hanced likelihood that r-process was successful in their PNS 
winds. 

We should caution that the comparison between our calcu- 
lations and the thresholds of HWQ in FigureQ~2]was made for 
Tdyn = 5 ms, which was the most rapid outflow HWQ consid- 
ered; this is, however, almost an order of magnitude longer 
than the dynamical time associated with some of our solu- 
tions. Constraints similar to HWQ at Td yn ~ 3 ms are obtained 
by Meyer & Brown (1997). A somewhat different thresh- 
old is obtained by Sasaqui et al. (2006), who emphasize a 
previously neglected light-element r-process seed production 

11 This conclusion disagrees with those of Nagataki & Kohri (2001), who 
also considered the effects of rotation and magnetic fields on PNS winds; 
however, these authors were unable to consider fields larger than ~ 10 11 G 
because of the complicated critical point topology they encountered in com- 
puting more highly-magnetized wind solutions. 



channel. In fact, recent work by Meyer (2002) shows that 
the threshold for r-process nucleosynthesis may be modified 
at very short Td yn compared to standard expectations (i.e., Y" 
must be less than 0.5). In order to accurately assess whether r- 
process nucleosynthesis in proto-magnetar winds will proceed 
to the third abundance peak and beyond, or whether the nu- 
cleosynthesis that results from these outflows may be used to 
constrain the birth rate of proto-magnetars, a survey of nucle- 
osynthesis calculations should be carried out at short Td yn . In 
fact, the critical question of whether this modest S a , very-low 
Td yn mode of r-process can reproduce the seemingly univer- 
sal solar abundance curve above Ba (e.g., Cowan et al. 2005) 
must ultimately be answered through detailed nucleosynthe- 
sis calculations on our wind solutions (such calculations are 
currently in progress) and by including a better treatment of 
Y e in the wind. The effects of GR not included in our calcula- 
tions will tend to increase S a and Y", on balance probably in- 
creasing the likelihood of successful r-process (Fuller & Qian 
1996; Cardall & Fuller 1997). 

T01 found that in NRNM winds S a oc (T dyn ) 2 ; thus, even 
though S" ex Lp 2 increases with time as the PNS cools, it is 
difficult for a NRNM wind that does not produce conditions 
favorable for third peak r-process (eq. [36)) at early times to 
enter the regime for a successful r-process at later times. We 
might expect a modification to the PNS (S a , Td yn ) evolutionary 
track in the case of a proto-magnetar because the dynamical 
time in a magnetically-driven PNS wind is no longer set solely 
by the neutrino heating. In order to address this question, 
we explored how (5 ,fl ) 3 /Tdy n changes as the PNS cools. As 
shown in Figure [TT] from our calculations at = 3.5 and 
B v = 2.5 x 10 14 G we find that (S a ) 3 /T dyn changes by less than 
a factor of 2 from those at L^^i = 8 and B v = 2.5 x 10 14 G; in 
particular, the peak value of (S a ) 3 / T dyn at P ~ 3 ms remains 
essentially unchanged. Thus the conditions for r-process in 
magnetically-driven PNS winds do not vary strongly with lu- 
minosity over the range we have explored. 

To conclude this section, we briefly consider what con- 
straints can be placed on r-process in proto-magnetar winds 
if they are to be capable of having a significant effect on the 
heavy element abundance evolution of the Galaxy. If most 
SNe produce r-process elements, the total r-process-rich wind 
mass ejected per SN must be ~ 10~ 6 - 10~ 5 M Q to account for 
the total Galactic yield (e.g., Qian 2000); hence, because the 
magnetar birthrate is ~ 10% of the total neutron star birthrate, 
at least ~ 10~ 5 - 10~ 4 M Q of r-process-rich material must 
be ejected per magnetar birth (this number is quite uncer- 
tain because the magnetar birthrate is uncertain; see Woods & 
Thompson 2004). Because r-process in proto-magnetar winds 
is only likely to be successful for P > 2 - 3 ms, rotation does 
not significantly enhance mass loss from the PNS (Fig. ITOb . 
As a result, the required yield of r-process rich material per 
magnetar birth is similar to our estimates for the total mass 
extracted in NRNM PNS winds (~ 10" 4 M Q ). Thus we con- 
clude that if proto-magnetar winds are the dominant source 
for Galactic r-process, the r-process probably must occur early 
in the PNS cooling evolution (in the first few seconds, even 
earlier than for normal PNSs) 12 . The former in part justifies 
our concentration on high L Pf in FigureQT]and means that the 
question of whether magnetar birth is a significant source for 

12 Of course, if the magnetar birthrate has been underestimated and is com- 
parable to the total NS birthrate then r-process could occur somewhat later in 
the PNS evolution. 
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Galactic r-process is especially sensitive to the early evolution 
of the magnetic field and radius of the PNS. 

It is also worth noting that on the basis of a comparison 
between the solar r-process abundance pattern with mete- 
oritic abundances of 129 I and 182 Hf, Qian, Vogel, & Wasser- 
burg (1998) (see also Wasserburg et al. 1996) suggest a di- 
versity of r-process sites. In particular, they argue for a site 
with high frequency (roughly the Galactic SN rate) that gen- 
erates the A w 195 nuclei and a site 10 times less frequent 
that produces nuclei near A w 135. In order to satisfy the 
observational constraints, the latter site must eject 10 times 
the mass of the high-frequency site, per event. Although we 
have not proven that magnetars produce an r-process in any 
mass range (we save a detailed investigation for a future work) 
the ~ 10% birth fraction of magnetars, the characteristically 
larger total ejected mass (Fig. IflOl ), and the very different 
thermodynamics of their winds relative to NRNM PNS winds 
make it tempting to associate proto-magnetar winds with the 
low-frequency enrichment events advocated by Qian, Vogel 
& Wasserburg (1998). 

4.3.1. Wave Heating 

QW show that by including a heating source in addition 
to neutrinos outside a few PNS radii, the wind entropy is in- 
creased, the dynamical time is reduced, and the chances for a 
successful r-process can be substantially improved. Such an 
extended heating mechanism operates above the sun, where 
convectively excited waves are believed to heat the extended 
solar corona (for recent work see Cranmer & van Ballegooijen 
2005). The neutrino cooling luminosity of a PNS also drives 
vigorous convection during the Kelvin-Helmholtz epoch (e.g., 
Burrows & Lattimer 1986; Burrows, Hayes, & Fryxell 1995). 
It is likely that a fraction of the convective energy flux will be 
deposited into outgoing waves, which will then propagate into 
the PNS atmosphere and deposit their energy and momentum 
on a length scale of order a few PNS radii. The relative im- 
portance of hydrodynamic and MHD wave excitation likely 
depends on the magnetic field strength of the PNS. Here we 
focus on heating by MHD waves in the magnetospheres of 
strongly magnetized PNSs (Suzuki & Nagataki 2005), though 
hydrodynamic wave excitation may be important as well (e.g., 
Burrows et al. 2006a,b). If the energy flux in MHD waves at 
the PNS surface is F W (R V ) = (1 /2)pv\(5B /B r ) 2 \ Rv , we esti- 
mate that the total wave heating Q w « AttR^F^R^ is given 
by 



10 45 (fii 5 )- ( - - ) ergs s 



(38) 



where 5B ll = 8B{R V ) is the amplitude of the waves excited at 
the PNS surface (with a density p(R„) ps 10 12 g cm" 3 ). The 
ratio of the wave heating in equation d38l to the total neutrino 
heating Q v = 4-n r 2 pq v dr in the absence of rotation can be 
approximated as 



^0.03(fi 15 ) 3 (^ 



P„51 



(39) 



where Q v as 4.4 x 10 49 (L Pti 5i/8) 2,8 ergs s" 1 is taken from our 
NRNM calculations and agrees reasonably well with the re- 
sults of T01 (their Table 1). Equation (l39l illustrates that, for 
efficient wave excitation, wave heating at early times may be- 
come important for field strengths > 10 14 - 10 15 G. Equation 



also appears to suggest that wave heating will become 
substantially more important as the PNS cools and L Pe drops; 
however, whether this in fact occurs is unclear because the en- 
ergy flux in waves and the surface amplitude will likely 
decrease with L„ r as the convective flux decreases. 

In order to quantify the effects that wave heating have on the 
entropy and dynamical time in rotating PNS winds, we con- 
sider a concrete model in which we add Alfven wave pressure 
and heating to our solutions. In the entropy equation (eq. J9)), 
this leads to an additional source term of the form 

2 



v A FwjRy) 
v r +v A xR v p 



— exp 



<r-R„) 

xRu 



(40) 



where the factor v A /(v r +VA) < 1 accounts for the work done 
by the Alfven waves. Equation ( f40b concentrates the total 
wave heating (Q w = An r 2 pq w dr) on a radial length scale 
w X^v- Radially propagating Alfven waves also exert a pres- 
sure on the fluid (P„), which contributes a term to the right 
hand side of the radial momentum equation (eq. 0) of the 
form (Lamers & Cassinelli 1999; Suzuki & Nagataki 2005) 

_ 1 dP v _ q w (SB) 2 3v r + v A 1 dp 
p dr 

where 



2(v r + v A ) 32-np v r + v A p dr 



—r, (41) 



(5B) = 



v r + v A 



(42) 



We consider a variety of dissipation lengths (R v x) an d surface 
wave amplitudes {8B U ) and assess the effects of Alfven wave 
heating on our wind solutions. For reasons discussed at the 
end of i j4.31 we are primarily interested in the effect that wave 
heating has on r-process at early times; thus, we concentrate 
on wave heating applied to high luminosity solutions. 

Our wave heating calculations are summarized in Table 
[3j Motivated by the r-process threshold of HWQ (eq. [36)), 
we quantify the improvement towards a successful r-process 
through an "improvement factor": 



Av(2w,x) ; 



[^0.5Mev/ 7 dyn]e w 
[^0.5Mev/ Td y n ]e w = 



(43) 



where So.5Mev is the entropy at the radius where T = 0.5 MeV. 
Figure QT] shows that for L Ppi51 = 8,B„ = 10 15 G, and P = 3 
ms (corresponding to the peak in Fig. fTTl i an improvement 
factor I w > 6 is required to exceed the HWQ threshold for 
third-peak r-process for Y° sa 0.4, while for a NRNM solu- 
tion at Lj7 r ,5i = 8 successful r-process requires 7 W > 25. We 
should note, however, that HWQ assumed adiabatic expan- 
sion through the a-process temperature range (T « 0.5-0.2 
MeV), while for large x wave heating is important at these 
radii and the expansion is not adiabatic. This may modify the 
7 W required for successful r-process. 

For P ~ 3(130) ms, B = 10 15 G, and L c „ 5 i = 8, we find that 
7 W > 6(25) requires an amplitude 5B V /B U > 0.5-0.7 and a 
dissipation length x ~ 3. 13 This wave amplitude corresponds 
to a total wave heating > 4 x 10 49 ergs s" 1 , which is com- 
parable to the total neutrino heating for this solution and is 

13 For very small x> the wave heating effectively acts as an increase to the 
neutrino luminosity, which decreases the asymptotic entropy S" (recall that 
- absent wave heating - the entropy decreases with increasing luminosity). 
For large \, on the other hand, the heating is concentrated outside the T = 0.5 
MeV radius and, while S" increases substantially, the entropy at T = 0.5 MeV 
remains relatively unaffected. 
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~ 10 of the neutrino luminosity at early times. The fact 
that a comparable wave energy is required in the slow and 
rapidly rotating limits appears somewhat coincidental. For the 
non-rotating solutions, the wave heating only reduces the dy- 
namical time to ~ 2 ms, never reaching the regime of Td yn < 
1 ms obtained in rotating PNS winds. The entropy increase 
at T = 0.5 MeV is, however, larger in the non-rotating case 
because the T = 0.5 MeV radius is at larger radii (and thus 
So.5MeV is closer to S"). Although the actual wave dissipation 
mechanism in the PNS atmosphere is uncertain, we note that 
for initial amplitudes of SB U /B V ~ 0.5, conservation of ac- 
tion implies that the waves are nonlinear with 5B/B r ~ 1 at 
~ 3 - 4R V , similar to the dissipation lengths (x) that are op- 
timal for r-process. Also note that for more rapidly rotating 
solutions (Pssl ms), significantly more wave heating is re- 
quired for successful r-process: even 8B V « B v is insufficient 
at B sa 10 15 G. This is because the total neutrino heating itself 
increases at P ~ 1 ms (p is larger in the neutrino heating re- 
gion due to magneto-centrifugal support), so a given amount 
of wave heating has less of an effect on the solution. 

For lower neutrino luminosities, smaller gw can lead to suc- 
cessful r-process. For example, at Lp 51 = 3.5 and B = 2.5 x 

10 14 G, we find that, for P = 3(130) ms, g w > 10 49 (5 x 10 48 ) 
ergs s" 1 and \ ~ 3 is required to eclipse the HWQ r-process 
threshold for Y" = 0.4. 14 The required wave heating in this 
case is again comparable to the total neutrino heating, but is 
only ~ 6 x 10~ 4 (3 x 10~ 4 ) of the neutrino luminosity of the 
PNS. At even lower L v , a yet smaller fraction of the neutrino 
luminosity in wave heating would be capable of yielding suc- 
cessful r-process. However, as discussed in the previous sec- 
tion, constraints on the r-process rich material required per 
magnetar-birth imply that the r-process must be successful at 
relatively high neutrino luminosities. We thus conclude that 
wave heating leads to successful r-process if^lO^-lO^of 
the PNS's neutrino luminosity emerges in wave power at early 
times in the Kelvin Helmholtz epoch, and that this required 
level of wave heating is essentially independent of magnetic 
field strength and rotation rate (for P > 2 ms). 

Although the mechanism and formula we describe here are 
appropriate to proto-magnetars with large magnetic fields, it is 
possible that hydrodynamic (as opposed to MHD) wave heat- 
ing is important and generic to normal PNS birth. In this case, 
waves may be generated by convective motions as the PNS 
cools, or via global modes of the PNS similar to those ob- 
served by Burrows et al. (2006a,b), that persist into the cool- 
ing epoch. As in the MHD case, a fraction of the total neutrino 
luminosity (10~ 4 — 10~ 3 ) must emerge in wave power over tkh 
to produce conditions suitable for the r-process. Because such 
a mechanism may operate generically (and not just in proto- 
magnetars) the requirement of an early-time r-process high- 
lighted above for the proto-magnetars and in §4.3 is somewhat 
alleviated. 

4.4. Additional Applications 

14 Suzuki & Nagataki (2005) found that for non-rotating solutions with 
Alfven wave heating, r-process was successful for Q w ~ 2 X 10 48 ergs s 
and x ~ 5 - 10 at Lp e $\ ~ 1. Our results are similar to theirs although their 
required wave heating is somewhat larger than an extrapolation of our results 
would suggest. However, they used different mean neutrino energies than we 
have assumed and obtain M ~ 2 X 10 -6 Mq s _1 , roughly 3 times greater than 
we predict from equation (Tj) at L B 5] = 1. This larger mass-loading may 
explain why their calculations required somewhat more wave heating to enter 
the regime of successful r-process. 



While most magnetar formation probably results from the 
core collapse of a massive star, magnetars may also form 
through the accretion-induced collapse (AIC) of a white dwarf 
(Nomoto et al. 1979; Usov 1992; Woosley & Baron 1992). 
Rapid rotation will automatically accompany AIC due to the 
accretion of mass and angular momentum prior to collapse, 
and a strong field may accompany the final stages of the 
PNS's contraction, amplified through either magnetic flux- 
freezing of the progenitor white dwarf's field or via dynamo 
action (Duncan & Thompson 1992). The properties of PNSs 
formed following AIC would thus be very similar to those 
of the proto-magnetars that we have considered in this pa- 
per and the resulting proto-magnetar wind would be accu- 
rately modeled using our calculations. AIC may thus give 
rise to a LGRB with properties similar to those considered in 
£|4.2| However, LGRBs from AIC will not produce significant 
amounts of Ni (Woosley & Baron 1992; Dessart et al. 2006) 
and hence will not be associated with a simultaneous hyper- 
energetic Type-Ic SN, as in, e.g. GRB 980425/SN1998bw or 
GRB 030329/2003dh (Galama et al. 1998; Hjorth et al. 2003; 
Stanek et al. 2003). Instead, AIC should be associated with 
a class of SN-less LGRBs like GRB 060505 and 060614 
(Fynbo et al. 2006; Gal- Yam et al. 2006; Delia Valle et 
al. 2006). A prediction of this model is that these LGRBs 
should be associated with both relatively old (few Gyr) and 
relatively young (^100 Myr) stellar populations. If some SN- 
less LGRBs are found to be associated with old stellar popu- 
lations, it would strongly support the AIC interpretation. 

Although we have focused on magnetized PNS evolution 
in this paper, the physical conditions we have considered are 
quite similar to those expected in a neutrino-cooled accre- 
tion disk surrounding a newly-formed black hole in the "col- 
lapsar" model for LGRBs (Woosley 1993; Paczynski 1998; 
MacFadyen & Woosley 1999). For concreteness, consider the 
properties of a thin accretion disk with an accretion rate M acc = 
0.2M© s" 1 and viscosity parameter a = 0.01 at a fiducial ra- 
dius of « 100 km, or approximately 10 Schwarzschild radii 
above a non-rotating black hole of mass 3M Q . According to 
the calculations of Chen & Beloborodov (2006), at these radii 
(where P w 10 ms) the disk is optically thick to neutrinos, 
the surface temperature of the disk is very similar to that of 
a PNS neutrinosphere (ss 2 — 3 MeV), and the total neutrino 
luminosity of the disk is L v ~ 0.04M acc c 2 » 1.4 X 10 52 ergs 
s" 1 (their Fig. 18), comparable to that of a PNS at early times 
(note that at lower accretion rates the disk will be optically 
thin to neutrinos). If this disk were threaded with a large-scale 
poloidal field with a strength corresponding to that expected 
for MRI turbulence, B w 10 14 - 10 15 G (plasma R w 10- 100), 
the physical conditions would indeed resemble those in proto- 
magnetar winds. Thus neutrino-magneto-centrifugal driving 
may be important in setting the mass-loading and energy loss 
rate in outflows from collapsar disks. Indeed, Levinson (2006) 
has calculated the mass loading of neutrino-driven outflows 
in general relativistic MHD for conditions anticipated in col- 
lapsar disks, finding qualitatively similar results to those dis- 
cussed in this paper for PNS winds. 

As a final context in which our calculations may be rele- 
vant, we note that a magnetized accretion flow like that con- 
sidered above for collapsars or a short-lived, rapidly rotating 
proto-magnetar may be formed following the merger of a NS- 
NS binary (Paczynski 1986; Eichler et al. 1989; Rosswog & 
Liebendorfer 2003). Our calculations would also describe 
magnetized outflow from these objects, but additional work 
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is needed to explore this application in more detail. 

5. CONCLUSIONS 

We have solved the one-dimensional non-relativistic 
neutrino-heated MHD wind problem in order to study the ef- 
fects that magnetic fields and rotation have on PNS wind evo- 
lution following the launch of the SN shock but prior to the 
end of the Kelvin-Helmholtz cooling epoch. Figure Q] sum- 
marizes the physical regimes of PNS winds in the presence of 
rotation and magnetic fields. We map our monopole solutions 
onto the axisymmetric, relativistic dipole calculations of B06, 
thus taking into account the effects that neutrino-driven mass 
loss have on the fraction of open magnetic flux and the PNS 
spin-down rate. Our primary conclusions are as follows: 

• We identify three types of PNSs based on the domi- 
nant character of their energy loss during the Kelvin- 
Helmholtz epoch (see Figure|9]): 

(1) For slow rotation and low magnetic field strengths 
(B* p < lO" 2 ^, where B di P = lO 14 ^ G is the 
surface dipole field strength of the PNS and 
P = l^ms ms is its birth period), a neutrino- 
heated, thermally-driven wind dominates energy 
and mass loss from the PNS; most radio pulsars 
were probably this type at birth. 

(2) For larger rotation rates and field strengths 
(1(T*J* < B f» < 2/> 2 s exp[2P" 2 ]), a non- 
relativistic, magnetically-driven wind dominates 
during the Kelvin-Helmholtz epoch. Most ob- 
served Galactic magnetars (~ 10% the birthrate 
of radio pulsars) were probably of this type at 
birth. For B ~ 10 15 G and P < 2 ms, greater 
than 10 51 ergs can be lost to a non-relativistic, 
magnetically-driven outflow during tkh, and, 
over a broad range of initial spin period, the 
energy extracted is many times larger than from 
a slowly rotating PNS (see Fig. [9). The outflow 
from this type of PNS is likely to be collimated by 
magnetic stresses and the asymmetric injection of 
energy may be sufficient to generate an anisotropy 
in the morphology of the SNe remnant. 

(3) For rapid rotation and field strengths some- 
what larger than those observed from Galactic 
magnetars (Z?^ p > 2P I 2 ls exp[2P~ 2 ]), a relativistic, 
magnetically-driven wind dominates energy loss 
during the Kelvin-Helmholtz epoch. Although the 
birthrate of PNSs of this type is probably small (if 
they are produced at all), they may be capable of 
producing hyper-energetic SNe and long-duration 
gamma-ray bursts. 

• For P ss 3(1) ms and Bf p > 10 15 G the total mass loss 
during the Kelvin-Helmholtz epoch is enhanced by a 
factor of w 2(10 2 ) relative to a non-rotating PNS (Fig. 

• For initial PNS spin periods of P w 1 ms and mag- 
netic field strengths of Bf p » 3 x 10 15 - 10 16 G, we find 
that > 10 52 ergs of rotational energy is extracted on a 
timescale of 10-40 seconds and that the magnetiza- 
tion a of the outflow is ^ 0.1 - 1000. The energy, lu- 
minosity, timescale, and mass-loading of the late-time 



outflow are all consistent with those required to ex- 
plain long duration gamma-ray bursts (assuming effi- 
cient dissipation of magnetic energy into kinetic energy 
at large radii; e.g., Drenkhahn & Spruit 2002). In ad- 
dition, outflows from such PNSs have the property that 
the PNS rotational energy is extracted with a roughly 
uniform distribution in log(er) over the timescale tkh- 
For these modest a winds, especially at early times, 
energy loss from the PNS is enhanced relative to pure 
force-free spindown because of addit ional open mag- 
netic flux (see the discussion in 34. It . Thus, a signif- 
icant portion of the PNS rotational energy can be ex- 
tracted in just a few seconds following the launch of 
the SN shock, perhaps sufficiently rapidly to increase 
the nucleosynthetic yield of the SN (e.g., 56 Ni). 

• Winds from PNSs with Bf v > 10 15 G and P w 2 - 10 
ms produce conditions almost an order-of-magnitude 
more favorable for third-peak r-process nucleosynthe- 
sis in the space of (S fl ) 3 / T dyn than do winds from more 
slowly-rotating, less-magnetized PNSs (see Fig. fTTb . 
For these rotation rates, the asymptotic entropy is sim- 
ilar to that of a non-rotating, non-magnetized wind, 
while the dynamical time is significantly reduced by 
magnetic acceleration (see Fig. [T2l . The very differ- 
ent thermodynamic properties of rapidly rotating proto- 
magnetar winds (relative to non-rotating PNS winds) 
may contribute to the inferred diversity of r-process 
sites (e.g., Qian et al. 1998). 

• Heating by outgoing hydrodynamic or MHD waves 
may be important in PNS winds. We find that if > 
10~ 4 - 10~ 3 of the PNS's neutrino luminosity emerges 
in wave power at early times in the Kelvin-Helmholtz 
epoch, then r-process is successful in PNS winds. This 
conclusion is relatively independent of magnetic field 
strength and rotation period for P > 2 ms. 
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PNS Wind Properties at Lp„ = 8 x 10 51 ergs s -1 , 
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TABLE 2 

PNS Wind Properties at B v = 2.5 x 10 14 G,R„ = 10 km, andM= 1.4 M 
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TABLE 3 

PNS Wind Properties with Wave Heating at L^.m = 8 and B v = 10 15 G 
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